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Chapter  1 


Introduction 


The  rotor  of  a  helicopter  can  be  modeled  in  different  ways.  For  actuator-disk  theories, 
there  is  a  lift  or  thrust  vector  perpendicular  to  the  plane  of  the  rotor,  pointing  in  the 
direction  of  Uft.  In  more  advanced  theories,  the  lift  is  taken  perpendicular  to  the  vortex 
sheet  that  trails  each  blade.  Simple  physics  state  that  the  reaction  to  this  thrust  is  an  induced 
velocity  in  the  opposite  direction.  The  induced  velocity  has  a  kinetic  energy  that  must  be 
provided  by  the  rotor.  This  energy  produces  the  induced  power  loss  of  a  rotor.  The  study 
of  induced  power  on  a  helicopter  rotor  has  been  developed  from  a  variety  of  perspectives 
along  the  years.  Induced  power  can  never  be  eliminated,  but  it  can  be  minimized.  An 
optimum  propeller  has  minimum  induced  power. 

The  study  of  induced  power  on  a  helicopter  can  be  divided  in  two  main  flight  regimes: 
axial  flow  and  skewed  flow.  The  first  regime  corresponds  to  the  case  for  which  the  flow 
coming  into  the  rotor  is  perpendicular  to  the  rotor  plane.  A  sketch  of  the  flow  and  rotor  can 
be  seen  in  Figure  1.2.  For  axial  flow  the  skew  angle,  y,  is  zero.  This  happens  when  the 
helicopter  is  in  hover  or  ascending  vertically  through  the  air.  When  the  helicopter  flies 
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forward  the  flow  coming  into  the  rotor  is  no  longer  perpendicular;  and  thus  the  skew  angle 
is  non-zero  and  a  skewed  flow  analysis  is  required. 


For  minimum  induced  power  purposes,  hover  and  forward  flight  are  the  both  important 
cases  to  consider.  If  the  induced  power  can  be  reduced  in  hover  the  helicopter  would  need 
less  power  in  order  to  lift  from  the  ground,  resulting  in  additional  power  that  could  be  used 
to  lift  more  payload  or  increase  the  range  of  the  helicopter.  Reducing  the  induced  power  in 
forward  flight  contributes  considerably  to  an  increase  in  cruise  speed  and  a  further  increase 
in  range. 

For  both  axial  and  skewed  flows,  we  can  model  the  rotor  in  one  of  two  ways:  as  an 
actuator  disk,  or  including  a  tilt  on  the  lift.  An  actuator  disk  is  a  way  to  model  the  rotor  as  a 
circular  region  of  zero  thickness  that  can  withstand  changes  in  pressure  across  it.  The  main 
assumption  with  an  actuator  disk  is  that  the  lift  vector  is  perpendicular  to  the  disk.  A  real 
rotor  has  some  tilt  on  the  lift  vector,  so  the  rotor  can  also  be  modeled  to  account  for  this 
condition.  If  the  lift  is  tilted,  the  wake  velocity  must  be  divided  into  two  components 
(normal  and  tangential  to  the  disk)  because  the  tilt  creates  a  swirl  velocity.  Some  changes  to 
existing  theories  need  to  be  made  to  account  for  the  tilt. 


Finally,  the  number  of  blades  considered  in  the  study  will  contribute  to  changes  in  the 
behavior.  Theories  must  be  developed  for  an  infinite  number  of  blades  (rotor  as  a  disk)  and 
for  a  finite  number  of  blades  (more  realistic).  It  is  of  importance  to  include  varying  numbers 
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of  blades  to  be  able  to  show  and  understand  how  varying  the  number  of  blades  affects  the 
induced  power. 


Table  1.1  shows  the  division  of  the  main  cases  that  will  be  studied.  Basically  we  have 
axial  or  skewed  flow,  depending  on  if  the  rotor  is  flying  through  air  vertically  or  on  forward 
flight  motion.  The  rotor  in  each  case  can  be  modeled  either  as  an  actuator  disk  with 
perpendicular  lift  vector  or  a  tilt  on  the  vector  can  be  included.  Lastly,  we  can  consider  a 
finite  or  an  infinite  number  of  blades. 


Previous  theories  predict  the  optimum  inflow  for  each  of  the  cases  in  axial  flow.  These 
theories  will  serve  as  references  to  determine  the  accuracy  of  the  present  model  and  better 
understand  the  changes  in  induced  power  for  the  different  cases  described  previously.  Once 
the  model  is  verified,  results  can  be  developed  for  skewed  flow  (or  forward  flight)  for 
different  configurations  with  the  security  that  the  approach  can  be  applied  for  such 
conditions. 


Table  1.1:  Cases  studied  using  finite  state  methods 


FLOW 

ACTUATOR  DISK 

TILTED  LIFT 

Infinite  Q 

Finite  Q 

Infinite  Q 

Finite  Q 

AXIAL 

Momentum  Theory 

Prandtl  Solution 

Betz  Distribution 

Goldstein  Distribution 

SKEWED 

Momentum  Theorj^ 

Finite-state 

Finite-state 

Finite-state  (theorj^  only) 

Several  alternative  approaches  can  be  used  to  model  the  downwash  in  the  plane  of  a 


rotor.  Techniques  such  as  CFD  and  free-wake  models  can  be  applied  successfully  for 
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different  applications  but  do  not  provide  best  results  for  real-time  simulations  due  to  the 
complexity  of  the  computations.  On  the  other  hand,  finite-state  models  provide  good 
results  for  preliminary  design  and  can  be  used  in  real-time  simulations.  Over  the  past  25 
years,  finite-state  inflow  equations  have  been  developed  to  describe  the  flow.  The  inflow  is 
expanded  in  terms  of  shape  functions,  resulting  in  a  matrix  set  of  ordinary  differential 
equations  that  describe  the  flow.  (See  Refs.  [l]-[9]).  The  development  of  these  models  has 
now  culminated  with  the  most  general  form  of  the  equations,  with  previous  versions  being 
special  cases  of  the  general  one. 

Although  the  rotor  blade  may  be  moving  at  transonic  speeds  near  the  tip,  the  rotor  wake 
is  generally  determined  by  incompressible  potential  flow  theory.  The  solution  of  the 
incompressible  potential  flow  equations  (and  of  their  derivative,  the  vorticity  transport 
equations)  has  generally  been  deemed  sufficient  to  define  the  wake  geometry  and  induced 
flow  of  rotors.  Thus,  most  models  are  based  on  these  equations,  including  finite-state 
models.  However,  previous  finite-state  applications  do  not  compute  the  flow  bellow  the 
plane  of  the  rotor  disk. 


One  purpose  of  this  dissertation  is  to  oudine  the  features  of  a  general  set  of  inflow 
equations,  including  their  derivation  and  special  cases  that  will  give  all  three  components  of 
velocity  everywhere  in  the  flow  field  including  within  the  wake.  Thus,  this  dissertation  can 
form  the  basis  for  any  version  of  the  model  for  anyone  interested  in  using  it.  The  present 
dissertation  also  provides  insights  into  the  limitations  and  accuracy  both  of  the  most  general 
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case  and  of  the  special  cases.  The  general  model  will  then  be  reduced  to  some  of  the  special 
cases  and  applied  to  obtain  the  induced  power  for  an  optimum  propeller. 


1.1  Previous  Work 


Many  attempts  have  been  made  with  the  purpose  of  modeling  and  understanding  the 
behavior  of  a  helicopter  induced  flow  in  order  to  reduce  the  negative  impact  on  the 
performance  of  the  rotor.  Previous  research  approaches  vary  in  complexity,  but  basically 
differ  from  each  other  depending  on  the  characteristics  of  the  induced  flow  distribution 
assumed  for  the  creation  of  the  theories.  Previously  developed  research  will  be  reviewed 
chronologically,  dividing  it  in  what  is  known  as  classical  work  and  more  recent 
developments. 

1.1.1  Classical  Work  on  Optimally  Efficient  Propellers 


In  1929,  Glauert  [15]  defined  the  minimum  induced  power  required  for  an  actuator  disk 
to  create  thrust  by  simple  work-energy  principles.  This  simple  Glauert  momentum  theory 
can  be  used  to  compute  the  minimum  possible  induced  power  for  a  given  flight  condition. 
In  addition,  he  included  the  geometry  of  the  flow  in  the  blade  coordinate  that  can  be  used  to 
obtain  a  relationship  between  induced  power  and  induced  velocity. 
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Originally,  momentum  theory  [16],  [17]  was  developed  to  describe  a  mathematical  model 
of  an  ideal  propeller  or  helicopter  rotor.  It  utilizes  the  basic  fluid  mechanics  conservation 
laws  to  predict  the  behavior  of  the  flowfield  of  a  rotor,  modeled  as  an  infinitely  thin  disk. 
The  simplest  version  of  the  rotor  as  a  disk  assumes  a  constant  velocity.  Therefore,  the 
induced  velocity  field  is  considered  to  be  uniform.  The  theory  provides  a  relationship 
between  power,  the  radius  of  the  rotor  and  the  induced  velocity.  This  theoretical  model  can 
be  applied  for  optimum  induced  power  considerations  for  the  case  of  axial  flow  for  an 
actuator  disk  with  an  infinite  number  of  blades.  The  results  provide  only  an  approximation, 
and  do  not  produce  results  that  can  be  used  to  design  the  rotor  blades,  since  the  use  of  a 
disk  considers  no  loss  of  thrust  at  the  blade  tips,  and  the  vortices  and  losses  of  each  blade 
are  therefore  not  accounted  for.  Due  to  the  limitations  of  he  uniform  momentum  theory,  a 
different  approach  to  the  same  problem  was  proposed  and  is  known  as  blade-element, 
momentum  theory. 

Blade-element  theory  [16],  [17]  assumes  that  each  element  of  the  rotor  is  independent  of 
the  other  elements  in  terms  of  lift  and  induced  flow.  The  total  thmst  and  the  torque  of  the 
rotor  are  obtained  by  integration  of  loads  along  the  span.  Because  the  blade  tip  does  not 
depend  on  its  location,  a  correction  factor  must  be  introduced.  Betz  and  Prandd  [18]  used 
combined  momentum  and  blade-element  theory  to  obtain  an  approximation  to  include  the 
effect  of  tip  losses.  Tip  losses  are  thrust  losses  due  to  the  fact  that  the  edge  of  vortex  sheet 
has  special  edge  effects.  Prandd  included  a  correction  to  account  for  a  finite  number  of 
blades,  thus  enabling  the  development  of  results  for  an  actuator  disk  in  axial  dow  and  a  finite 
number  of  blades. 
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The  optimum  induced  flow  distribution  for  an  actuator  disk  with  an  infinite  number  of 
blades  is  uniform  flow.  However,  when  the  lift  vector  is  taken  as  perpendicular  to  the 
vortex  sheets  (rather  than  perpendicular  to  the  disk),  the  ideal  power  is  no  longer  obtained 
by  uniform  flow.  Betz  determined  that  the  minimum  power  is  obtained  when  the  induced 
flow  at  the  individual  blades  is  such  that  the  vortex  sheets  remain  along  helical  paths.  For  an 
infinite  number  of  blades,  it  follows  that  the  pressure  field  must  follow  this  same  shape. 

The  effect  of  a  finite  number  of  blades  incurs  an  additional  loss  in  wake  energy  due  to  the 
individual  vortex  sheets  from  each  blade.  In  1929,  Goldstein  [19]  worked  out  the  exact 
effect  for  optimized  rotors.  Prandtl  developed  a  correction  factor  that  agrees  very  well  with 
Goldstein’s  exact  solution  for  moderate  cUmb  rates. 


1.1.2  Classical  Work  on  Inflow  Theories 


In  1981,  Pitt  and  Peters  [1]  developed  a  set  of  differential  equations  to  obtain  the 
dynamic  inflow  derivatives  for  a  helicopter  rotor  with  an  unsteady  loading  and  induced  flow 
distribution.  The  theory  relates  rotor  loads  and  the  transient  response  of  the  rotor  induced 
flow  field.  The  dynamic  inflow  theory  was  developed  to  obtain  induced  inflow  solutions 
that  are  more  realistic  than  those  of  momentum  theory  and  which  provide  an  expression  of 
the  induced  velocity  distribution  for  a  helicopter  rotor.  The  Pitt-Peters  model  uses  three 
terms  to  describe  the  inflow  distribution.  The  model  is  in  use  in  a  wide  variety  of  companies 


and  institutions  that  use  it  for  stability  and  handling  quality  applications.  There  are  however 
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some  limitations  to  the  theory,  as  it  gives  only  a  very  low-order  approximation  to  the 
induced  flow  fleld  of  a  helicopter  rotor. 

He  and  Peters  [4],  [5],  [6]  further  developed  the  theory  to  obtain  a  generalized  dynamic 
wake  model  that  included  higher-order  terms.  The  unsteady  induced-flow  theory  can  be 
appUed  in  stability,  vibration,  control  and  aeroelastic  studies.  The  theory  is  based  on  an 
acceleration  potential  for  an  actuator  disk.  The  Pitt-Peters  model  included  only  the  zeroeth 
and  first  harmonics  with  one  radial  shape  function  for  each  harmonic.  The  He  theory 
expands  this  to  include  a  higher  number  of  harmonics  and  of  radial  shape  functions  for  each 
harmonic.  The  pressure  and  the  induced  flow  are  expressed  in  a  polynomial  distribution 
(proportional  to  Legendre  functions)  radially  and  in  terms  of  a  Fourier  series  azimuthaUy. 
By  the  use  of  higher  harmonics,  the  flow  can  be  studied  in  detail  for  other  applications 
instead  of  being  limited  to  obtaining  the  average  flow  in  only  for  thmst,  pitch  moment  and 
roll  moment,  which  was  the  case  with  the  Pitt-Peters  model.  However,  there  are  some 
limitations  to  the  use  of  the  generalized  dynamic  wake  model  proposed  by  He  and  Peters. 
One  of  the  limitations  is  that  the  model  does  not  account  for  flit  on  the  lift  vector.  It  is  an 
actuator  disk  theory.  Also,  the  flow  in  the  He  model  is  not  deflned  in  all  three  directions  or 
everywhere  in  the  rotor  flow  field,  therefore  more  work  needed  to  be  done  to  be  able  to 
analyze  situations  were  the  flow  has  to  be  deflned  everywhere  and  the  use  of  all  components 
in  needed. 

Morillo  and  Peters  [7]  redefined  a  new  dynamic  inflow  model  to  overcome  the  kmitaflons 
of  the  Peters-He  model.  The  equations  that  define  the  flow  are  re-derived  from  first 
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principles  to  compute  the  whole  flow  fleld  with  all  of  the  components  for  the  volume  above 
the  rotor  disk.  The  model  uses  continuity  and  momentum  equations  (partial  differential 
equations)  in  three  dimensions  around  the  actuator  disk  to  model  the  flow.  The 
conservation  equations  are  then  transformed  to  a  set  of  ordinary  differential  equations  in 
state-space  form  by  the  application  of  Galerkin  method.  The  set  of  ordinary  differential 
equations  gives  a  relation  between  the  velocity  coefficients  and  the  pressure  coefficients 
(similar  to  the  one  for  the  Peters-He  model)  but  includes  an  expansion  in  the  influence 
coefficient  and  damping  matrices  that  allowed  for  the  three  components  of  velocity  (axial, 
radial  and  azimuthal)  to  be  determined  everywhere  on  the  upper  hemisphere  of  the  disk. 
Because  of  the  fact  that  the  method  has  the  same  matrix  stmcture  than  the  one  for  the  Pitt- 
Peters  and  Peters-He  models,  the  addition  of  the  improved  features  would  be  easily  applied 
to  codes  that  already  utilize  these  previous  versions  of  the  dynamic  wake  model.  However, 
this  theory  is  not  able  to  deflne  the  velocity  fleld  below  the  plane  of  the  rotor  disk.  See 
Figure  3.2. 

Makinen  and  Peters  [10],  [11]  expanded  the  Peters-He  as  well  as  the  Morillo-Peters  model 
to  account  for  lift  tilt  by  the  introduction  of  a  swirl  correction.  Different  swirl  corrections 
were  analyzed  and  used  to  perform  comparisons  with  the  results  provided  by  Goldstein. 
Makinen  was  able  to  obtain  accurate  results  for  a  flnite  number  of  blades  by  the  addition  of 
the  extra  kinetic  energy.  The  comparison  of  the  velocity  fleld  and  the  circulation  around  the 
propeller  blades  was  goof  over  most  of  the  blade  span.  He’s  model  provides  the  necessary 
fidelity  to  be  able  to  obtain  the  induced  power  for  a  helicopter  rotor  modeled  as  an  actuator 
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disk,  and  the  addition  provided  by  Makinen  allows  to  extend  the  theory  for  the  tilted  lift  case 
of  a  rotor  with  a  finite  number  of  blades. 

Different  inflow  theories  have  been  developed,  Rosen  [24]  proposed  an  approximate 
actuator  disk  model  of  a  rotor  in  hover  or  axial  flow  based  on  potential  flow  equations.  The 
model  is  derived  from  the  same  fundamental  equations  as  the  ones  used  in  the  Refs.  [1]  and 
[4]  model,  and  in  fact  his  theory  shows  the  solution  for  the  Pitt-Peters  model  as  a  special 
case.  The  assumptions  made  enable  the  model  to  include  wake  distortion  effects  in  addition 
to  the  original  dynamic-inflow  model,  but  as  the  models  shown  in  references  [1]  through  [8], 
the  solution  is  only  applicable  to  actuator  disk  cases. 


1.1.3  Motivation  of  Optimization  Work 


Harris  [20]  examined  a  large  set  of  data  for  rotor  power  coefficient  versus  advance  ratio. 
Figure  1.1  shows  the  thrust  for  a  helicopter  in  forward  flight.  If  the  profile  and  parasite 
power  are  subtracted  from  the  total  helicopter  power  (the  two  higher-ordinate  curves),  one 
is  left  with  the  net  induced  power  of  a  helicopter.  These  results  show  that,  for  moderate 
advance  ratios,  the  difference  between  the  measured  rotor  induced  power  and  the  theoretical 
prediction  for  minimum  induced  power  proposed  by  Glauert  are  considerably  large.  This 
discrepancy  implies  that  lifting  rotors  in  forward  flight  have  efficiencies  far  lower  than  the 
optimum  efficiencies  predicted  by  theory. 
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Figure  1.1:  Harris’  power  coefficient  with  advance  ratio  plot  for  a  helicopter. 


Following  Harris’  results,  Ormiston  [21],  [22]  performed  extensive  mns  with  RCAS 
(Rotorcraft  Comprehensive  Analysis  System)  to  try  to  determine  why  the  actual  results  were 
differing  from  the  ideal  results.  In  these  studies,  the  profile  drag  of  the  airfoil  was  assumed 
to  be  zero  so  that  the  induced  power  could  be  separated  out.  The  results  of  those  studies 
similarly  showed  a  large  difference  between  ideal  induced  power  and  the  actual  induced 
power  of  rotorcraft,  the  actual  being  several  times  larger  than  ideal. 


In  the  hopes  of  expanding  on  the  understanding  of  optimum  rotors,  other  research 
efforts  have  used  other  methods  to  determine  why  there  is  a  considerable  difference 
between  the  theoretical  predicted  minimum  induced  power  and  the  real  results  obtained 
from  different  sources.  The  use  of  computer  codes  such  as  EHPIC/HERO  (Evaluation 
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Helicopter  Performance  using  Influence  Coefficients/  Helicopter  Rotor  Optimization), 
which  provides  optimization  techniques  for  the  design  of  rotor  blades,  or  the  software 
CHARM  (Comprehensive  Hierarchical  Aeromechanics  Rotorcraft  Model)  [23]  — that 
models  aerodynamics  and  dynamic  of  rotorcraft  in  general  flight  conditions —  have  been 
applied  to  obtain  induced  power  for  a  variety  of  flight  conditions.  AU  results  show  the  same 
apparent  shortfall. 


1.2  Optimization  Problem  Statement 


Work-Energy  principles  indicate  that  the  induced  power  Pj  generated  for  a  lifting  rotor 
(i.e.,  the  power  that  does  not  perform  useful  work)  can  be  found  by  computing  the  shaft 
power  and  then  subtracting  the  work  done  on  the  vehicle 

Pj  =Ps-TV-HU  (2.1) 

where  T  is  the  thrust  perpendicular  to  the  disk,  V  is  the  rotor  velocity  in  the  T  direction,  H  is 
the  rotor  force  in  the  inplane  direction,  and  U  is  the  rotor  velocity  in  the  H  direction  (See 
Figure  1.2).  By  necessity,  the  magnitude  of  this  power  must  equal  the  power  that  is 
expended  into  the  kinetic  energy  of  the  induced  flow.  It  follows  that  simple,  Glauert 
momentum  theory  can  be  used  to  compute  the  minimum  possible  induced  power  for  a 
given  flight  condition. 
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Several  potential  sources  of  decreased  efficiency  can  be  identified  in  terms  of  the  physics 
of  an  actual  rotor  as  compared  to  an  ideal  actuator  disk.  First,  an  ideal  disk  produces  thmst 
perpendicular  to  the  disk  whereas  a  true  rotor  produces  a  tilted  thrust  vector  that  results  in 
swirl  velocity.  Therefore,  there  is  lost  energy.  Second,  an  ideal  rotor  has  an  infinite  number 
of  blades  whereas  a  true  rotor  has  a  finite  number.  The  fact  that  there  are  vortex  sheets 
coming  off  the  individual  blades  implies  an  upwash  outside  of  the  slipstream  that  further 
translates  into  lost  energy.  Third,  an  ideal  disk  can  generate  an  arbitrary  lift  distribution  over 
the  span  and  azimuth.  An  actual  blade,  on  the  other  hand,  can  only  produce  lift  under  the 
constraints  of  both  allowable  blade  pitch  changes  and  of  the  limits  on  airfoil  lift  coefficients 
at  high  angles  of  attack.  One  goal  of  the  present  research  effort  is  to  determine  which  of 
these  contribute  to  the  drastic  increase  in  induced  power  and,  consequently,  what  changes  in 
rotor  hardware  (if  any)  might  address  the  issue. 


Figure  1.2:  Basic  illustrative  problem  for  a  rotor  [22]. 
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1.3  Present  Approach 


Finite-state  inflow  models  have  been  suggested  as  a  theoretical  basis  whereby  to  study  the 
reasons  for  this  efficiency  deflcit.  In  this  dissertation,  a  flnite-state  inflow  model  is  utilized 
to  formulate  the  optimum  circulation  and  inflow  distribution  for  rotors  in  axial  flow.  The 
results  show  that  a  formal  optimization  with  flnite-state  models  can  be  done  in  closed  form 
and  that  such  an  optimization  recovers  the  classical  uniform-flow  condition  (for  an  actuator 
disk  with  an  infinite  number  of  blades),  the  Prandtl  solution  (for  an  actuator  disk  with  a 
flnite  number  of  blades),  the  Betz  distribution  (for  a  lifting  rotor  with  an  infinite  number  of 
blades)  and  the  Goldstein  solution  (for  a  lifting  rotor  with  a  finite  number  of  blades).  Thus, 
it  should  be  possible  to  use  flnite-state  models  to  investigate  optimum  rotor  performance  in 
forward  flight. 


In  classical  momentum  theory,  the  optimum  performance  of  an  actuator  disk  is  well 
known  to  be  a  uniform  pressure  distribution  with  a  resulting  uniform  inflow.  However, 
propellers  and  rotors  are  not  actuator  disks.  In  particular,  the  lift  vector  for  rotors  is 
perpendicular  to  the  vortex  sheets,  not  to  the  disk.  As  a  consequence,  the  optimum 
distribution  for  pressure  and  inflow  for  a  true  rotor  is  the  Betz  distribution  rather  than 
uniform.  Furthermore,  real  rotors  have  a  flnite  number  of  blades  which  further  modifies  the 
optimum  distribution  to  be  that  of  Prandtl  for  an  actuator  disk  and  that  of  Goldstein  for  a 
rotor  with  a  flnite  number  of  blades. 
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Interestingly,  finite-state  inflow  models  offer  a  natural  framework  in  which  to  study 
optimum  inflow  distributions  because  —  in  these  models  —  both  inflow  and  thrust  are 
expanded  in  a  series  of  orthogonal  functions  that  make  the  optimization  problem  a  closed- 
form,  classic  case  of  minimizing  a  quadratic  cost  function  with  linear  constraint.  In  addition, 
recent  improvements  to  dynamic  wake  models  allow  for  both  the  tilt  of  the  thmst  vector 
and  the  resultant  swirl  velocity.  This  allows  finite-state  methodologies  to  be  rigorous  in  the 
context  of  optimization  of  induced  power. 

In  this  dissertation,  we  show  how  a  closed-form  optimization  of  induced  power  with 
finite-state  models  reproduces  the  momentum  result  for  an  actuator  disk,  the  Prandtl  result 
for  a  disk  with  a  finite  number  of  blades,  the  Betz  result  for  an  ideal  rotor,  and  the  Goldstein 
result  for  a  rotor  with  finite  number  of  blades.  This  success  gives  reason  to  believe  that  the 
methodology  can  be  used  in  the  future  for  applications  to  optimization  in  forward  flight. 

One  goal  of  this  research  is  to  derive  the  theoretical  background  to  obtain  results  for  an 
optimum  propeller.  The  model  will  provide  the  necessary  solutions  to  be  used  in  real 
applications  by  the  addition  of  the  resulting  expressions  to  the  general  code  already  in  use 
for  many  applications.  Results  are  shown  to  validate  the  theory  and  show  the  induced  power 
for  the  cases  where  theory  was  never  previously  developed.  The  optimization  provides  a 
general  solution  of  the  pressure  coefficients  for  a  general  rotor  blade  configuration. 
However,  it  can  be  applied  for  specific  cases  (specific  blades  or  configurations)  with  changes 
to  the  constraints  and  parameters  of  the  optimization. 


15 


Chapter  2 


Complete  Inflow  Theory  Formulation 


Finite-state  inflow  theories  have  been  used  in  order  to  define  the  wake  geometry  and  the 
induced  flow  of  rotors.  The  Pitt-Peters  model  and  the  higher  order  Peters-He  model 
represent  the  downwash  in  the  plane  of  the  rotor  being  expanded  in  terms  of  a  set  of  radial 
an  azimuthal  functions  that  describe  the  flow.  The  inflow  is  therefore  described  by  a  finite 
set  of  state  equations  that  are  driven  by  the  loading  on  the  rotor  blades. 

Up  to  this  point,  previous  developments  have  been  able  to  achieve  the  full  description  of 
the  wake  in  all  ranges  of  flight  regime  —  for  all  skew  angles  including  axial  and  edgewise 
conditions  for  all  three  components  of  the  flow  —  axial,  radial  and  azimuthal  — ,  and  in  all 
regions  around  an  actuator  disk  except  for  the  lower  hemisphere.  Figure  2.1  shows  the 
different  areas  around  the  actuator  disk,  where  regions  B  and  C  lay  on  the  lower  hemisphere. 
The  purpose  of  this  chapter  is  to  provide  a  set  of  ordinary  differential  equations  in  matrix 
form  that  describe  the  flow  for  all  ranges,  including  all  three  components,  and  for  all  regions. 
This  formulation  will  give  the  most  general  form  of  the  equations,  with  previous  versions 
being  special  cases  of  the  general  one.  It  is  thus,  a  complete,  nonlinear  induced  flow  theory 
for  rotors  in  incompressible  flow. 
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A 


Figure  2.1:  Regions  on,  off  and  around  an  Actuator  Disk 


2.1  Potential  Flow  Equations 


Conservation  of  mass  and  momentum  (partial  differential  equations)  in  three  dimensions 
around  an  actuator  disk  can  be  reduced  to  the  nonlinear,  potential  flow  equations  for 
incompressible  flow.  Ref.  [25]. 

V-V  =  0  (2.1) 

^  +  {v.v)v  =  -^  (2.2) 

0/  ^  P 

In  past  studies,  these  have  been  treated  in  various  nondimensional  schemes.  Here,  for 
clarity,  we  keep  everything  dimensional. 
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In  the  following  derivation,  we  will  first  lineari2e  about  a  steady  free-stream  velocity. 
Later,  we  will  come  back  and  consider  how  some  limited  nonknearities  can  be  introduced 
into  the  equations.  The  coordinate  system  in  Figure  2.2  shows  the  rotor  disk  with  x  in  the 
nominally  forward  direction,  y  into  the  paper  (the  advancing  side  of  the  rotor),  and  z  positive 
down. 


The  steady  free-stream  is  considered  to  be  of  magnitude  Voo  and  directed  opposite  to  a 
unit  vector  ^  that  is  skewed  with  respect  to  the  negative-^  axis  by  an  angle  positive 
towards  the  x  axis.  Thus,  we  have  for  the  free-stream  velocity, 

V  =  -Vj  +  v  (2.3) 

I  =  sin{z)i  -  cos{z)j  (2.4) 
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where  i,j,  k  are  unit  vectors  in  the  x,  j,  and  ^  directions,  respectively.  When  these 


definitions  are  put  into  the  momentum  equations  and  linearized,  the  resultant  equations 
become: 


V  •?  =  0 

The  boundary  condition  is  that  all  components  of  vmust  be  zero  far  upstream. 


(2.5) 

(2.6) 


2.2  Pressure  Potentials 


The  next  step  in  the  derivation  is  to  expand  the  pressure  of  the  linearized  equations  in 
terms  of  potential  functions.  Discontinuities  in  pressure  (or  in  its  derivatives)  can  only  occur 
on  the  disk.  Therefore,  we  utilize  the  solution  to  Laplace’s  equations  in  ellipsoidal 
coordinates,  defined  on  Appendix  B,  to  represent  the  pressures  in  the  flow  field.  This  then 
allows  the  introduction  of  time-varying  pressure  drops  or  mass  sources  at  arbitrary  locations 
on  the  rotor  disk.  The  solutions  that  we  include  are  those  that  die  out  at  infinity.  They  are; 

O:  =  P;  (v)^:  {iri)co,myf  (2.7) 


=  K  {y)Q” 


(2.8) 


where  V,  rj,  and  if/  are  the  ellipsoidal  coordinates  defined  in  Appendix  A  and  m  takes  the 
values  {0,  1,  2,  3,  4,  ...  }and  n  takes  on  the  values  {m,  m+1 ,  m+2,  m+3,  ...}. 
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It  is  interesting  to  see  how  these  potential  functions  behave  across  the  ^—0  plane  (where 
V  jumps  from  positive  to  negative).  In  the  plane  of  the  rotor  disk  but  off  of  the  rotor  disk, 
we  have  v=0.  Therefore,  the  functions  with  m+n  odd  (which  are  odd  functions  in  v)  are 
zero  on  that  plane;  and  the  functions  with  m+n  even  (which  are  even  functions  in  v)  have  a 
zero  normal  derivative  on  that  plane.  On  the  upper  surface  of  the  disk,  we  have  0<v<1\ 
and,  on  the  lower  surface,  we  have  -1  <v<0.  Thus,  V  jumps  across  the  disk.  Therefore,  the 
m=n  odd  potentials  have  a  jump  across  the  disk  (simulating  a  pressure  jump);  and  m+n  even 
functions  have  a  jump  in  normal  derivative  (simulating  a  mass  source). 

The  pressure  may  now  be  expanded  in  terms  of  these  potentials. 


P 

P 


00  00 

_y  y  t  ) 

\  n  n  n  n  J 

w=0 n^m 


(2.9) 


where  the  O’s  are  nondimensional  so  that  T  will  have  the  units  of  m^/sec^.  Note  that,  for 
m-^n  odd,  the  pressure  will  be  discontinuous  across  the  disk  since  (k)  is  odd  in  that  case. 
Similarly,  for  m-^rn  even,  the  normal  gradient  of  pressure  will  have  a  discontinuity  across  the 
disk.  Thus,  m-^rn  odd  corresponds  to  lift  on  the  disk,  and  m-^rn  even  corresponds  to  a  mass 
source  on  the  disk.  It  follows  that  the  following  relationships  connect  the  t’s  to  either  the 
pressure  drop  or  the  injected  mass  on  the  plane  of  the  disk: 


00  00 


AP  =  P,  -P,,  =2py  y  + 

lower  upper  r  /  j  /  j  \  ^  ^  ^  J  ri={) 


(2.10) 


m=0  «=«+!, z^+3,. 


ip  +P  ) 

.  ,  .  \  lower  upper  }  _ 

=  pAv  =  =  -2 


f  P  ^  «1  00 

J  m=0  n=m+2,m+A^... 


Zy 

\  n  n  n  n  J 


77=0 


(2.11) 


where  AP  is  the  net  pressure  in  the  upward  direction  across  the  disk  (negative  and  Am  is 


the  net  mass  per  unit  area  being  injected  into  the  flow  field  per  unit  time  at  the  disk 
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(implying  a  positive  ^  velocity  below  the  disk  and  negative  ^  velocity  above  the  disk.  Note 


that  only  m-^rn  odd  will  contribute  to  AP,  and  only  m+n  even  will  contribute  to  Am .  Since 
the  potential  functions  are  a  complete  set  of  orthogonal  functions  on  the  disk  for  either  mPn 
odd  or  mPn  even,  knowledge  of  the  pressure  drop  and  mass  injection  at  the  disk  at  any  point 
in  time  completely  determines  the  x’s  at  that  instant  in  time. 

2.3  Velocity  Potentials 


The  first  velocity  potentials  that  we  define  will  be  called  the  prime  potentials,  W.  These 
are  defined  as  the  integrals  along  a  streamline  of  the  pressure  potentials. 


(2.12) 


where  the  rotor  radius  R  is  added  to  keep  the  velocity  potentials  nondimensional.  The  above 
relationship  holds  both  for  the  cosine  terms  (me)  and  for  the  sine  terms  (ms).  The  flow 
velocities  are  then  represented  by  an  expansion  summation  in  the  gradient  of  these  velocity 
potentials. 

00  00 

»=RZZ("T'5'r+*T'5'.")  (2-«) 

m=0 n=m 


Note  that,  since  the  W  are  nondimensional,  since  the  R  is  added  for  scaling;  and,  since  the 
Laplacean  operator  has  units  of  one  over  length,  the  a”  and  h"'  will  have  units  of  velocity, 
m/ sec. 
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It  is  straightforward  to  show  that  these  expansion  gradients,  V'F ,  are  exacdy  the  velocity 


due  to  each  individual  pressure  potential  (P  in  the  static  case. 


dv  _  VP 


(2.14) 


Substitute  Equation  (2.9)  and  (2.13)  into  Equation  (2.6).  This  leads  to: 


00  00 


K.££(«:v®:'+vv<i>r)=|vx  £  (rr®r+c'®:')»'f  (2.15) 


m=0 n=m 


»=0 


Then,  based  in  the  definition  in  Equation  (2.12), 


T  T 

_  n  .  ym  _  n 

^  n  ^  n 


(2.16) 


Thus,  this  particular  choice  of  velocity  potential  will  give  exact  results  at  zero  frequency  even 
for  truncated  matrices. 


The  above  definition  of  velocity  potentials  involves  three  areas  of  the  flow  field  that  will 
be  important  to  identify.  First,  there  is  the  area  above  the  plane  of  the  rotor  disk  {^>0,  ^<0, 
0<v<1).  This  is  the  upper  hemisphere  where  one  need  only  be  concerned  with  positive 
values  of  V.  A  second  area  is  the  flow  within  the  wake  (i.e.,  the  flow  below  the  plane  of  the 
rotor  disk  along  streamlines  that  pass  through  the  rotor  disk,  region  C  in  Figure  2.1).  This 
area  defines  a  skewed  cylinder  that  encompasses  the  rotor  wake.  When  one  follows  such  a 
streamline  into  the  wake,  vis  discontinuous  across  the  rotor  disk  resulting  in  discontinuities 
either  in  the  potential  functions  or  in  their  derivatives.  This  invalidates  the  assumptions  of 
Equatons.  (2.12-2.13)  since  the  Laplacean  operator  cannot  be  brought  out  from  under  the 
integral  sign.  Stricdy  speaking,  the  flow  is  not  potential  flow  within  the  wake,  giving  issues 
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in  finding  the  induced  flow  there.  The  third  area  of  the  flow  to  be  defined  is  the  flow  below 
the  plane  of  the  rotor — but  along  streamlines  that  do  not  intersect  the  rotor  disk.  The 
velocity  potentials  are  valid  in  that  area,  which  is  a  region  governed  by  potential  flow  (i.e., 
region  B  in  Figure  2.1). 


Unfortunately,  there  is  presendy  no  known,  closed-form  expression  for  these  velocity 
potentials  in  Equation  (2.12)  at  arbitrary  wake  skew  angles.  Thus,  in  the  above  formulation, 
some  preprocessing  must  be  done  to  compute  the  gradients  of  IP  at  points  in  the  flow  field 
for  which  the  velocity  is  desired.  In  order  to  overcome  this  limitation,  another  set  of 
velocity  potentials  is  defined,  the  derived  potentials,  y/  .  They  reduce  to  the  prime  potentials 
•P  for  the  special  case  in  which  X~0-  order  to  find  these  functions,  one  must  find 
potential  functions  such  that 


jvn*  jxu'” 


(2.17) 


Reference  [6]  was  able  to  find  these  functions  in  the  following  form: 


^  ^»+l  '  hn 


(2.18) 


where 


cr.  = 


P;  /(2»  +  l)(2»  +  3)  (»  +  l)"-^" 


m 


(2.19) 


and 
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k:  = 


H”  = 


h: 


\2J 

(»  +  /^  +  l)!!(»  — — l)!! 
(»  +  —  m)U 


(2.20) 


which  are  valid  for  the  upper-half  plane  (^<0)  with  the  exception  of  the  case  /^=».  For  that 
case,  i^is  infinite  and^f/T^J  and  PfvJ  are  undefined  for  the  subscript  less  than  the  superscript 


Reference  [9]  completed  the  formulation  for  these  special  cases  so  that  now  all  of  these 
derived  potential  functions  are  known  for  all  superscripts  and  subscripts.  For  the  case 
but  not  equal  to  zero,  the  derived  potentials  become: 

'i'r =[v'c, (•')£:;, ('>?)+f;c, (•')£:-, (2.21) 

'j'r (2.22) 

where  P*.  and  Q”  .  are  defined  as 

m—\ 
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(2.24) 


For  the  special  case  m—n—0,  the  velocity  potential  is: 

f  1  4 


n 


\-r]tavL 


1 

\T1) 


2  1  2 

- ln|l  +  v| - ln|l  +  ;7"|  +  -ln|Z^^  I  v  >0  (2.25) 

71  71  ''  71 
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where  is  a  large  number  representing  the  radius  to  which  the  ^  integral  is  taken.  In 
reality,  the  integral  is  infinite  when  taken  to  infinity.  However,  it  is  only  the  gradients  of  the 
potential  functions  that  are  required  — not  the  potential  functions  themselves —  and  one  can 
always  add  a  constant  to  a  potential  function.  Therefore,  the  above  definition  is  tractable  for 
practical  solutions.  For  the  derived  velocity  potentials  below  the  plane,  one  can  consider 
that  the  velocity  potentials  with  m-^rn  odd  must  be  symmetric  with  respect  to  ^=0;  while 
functions  with  m-^rn  even  must  be  asymmetric  to  that  plane.  Equation  (2.18).  That  defines 
the  potentials  in  the  lower  hemisphere  outside  of  the  wake,  region  B.  As  with  the  prime 
potentials,  potentials  within  the  wake  will  require  special  attention  and  will  be  considered 
later  in  this  chapter. 

Since  the  velocity  potentials  all  satisfy  Laplace’s  equation,  it  follows  that  the  velocities, 
themselves,  will  satisfy  continuity.  Equation  (2.1).  This  leaves  only  the  momentum  equation 
as  the  governing  equation  for  the  velocity  expansion  coefficients.  Therefore,  the  problem  of 
obtaining  a  finite-state  wake  model  is  reduced  to  the  problem  of  representing  the 
momentum  equations  in  finite-state  form. 

2.4  Galerkin  Approach 


The  next  step  in  the  derivation  is  to  use  the  Galerkin  approach  (described  on  Appendix 
C)  to  transform  the  momentum  equation  into  a  set  of  ordinary  differential  equations  for  the 
expansion  coefficients.  To  this  end,  the  pressure  expansions.  Equation  (2.9),  and  the  prime 
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velocity  potential  expansions.  Equation  (2.13)  are  substituted  into  Equation  (2.6).  The  test 
functions  are  taken  to  be  the  gradient  of  the  pressure  potentials.  These  test  functions  are 
dotted  with  the  momentum  equation  and  integrated  over  the  upper  hemisphere  (above  the 
plane  of  the  rotor  disk)  to  obtain  equations  of  motion.  (This  avoids  the  problems 
mentioned  earlier  about  discontinuities  across  the  rotor  disk  and  in  the  wake.)  It  should  also 
be  noted  that,  since  the  upper  hemisphere  is  chosen,  there  are  no  velocity  boundary 
conditions  across  vortex  sheets,  due  to  the  fact  that  discontinuities  only  occur  in  the  wake, 
which  is  below  the  plane  of  the  disk.  Thus,  either  the  prime  potentials  or  the  derived 
potentials  are  appropriate  for  a  Galerkin  approach.  As  a  result  of  integrating  over  the  upper 
hemisphere,  however,  the  velocity  solution  will  only  be  valid  in  the  upper  hemisphere. 
Later,  the  velocity  field  below  the  rotor  plane  will  be  discussed. 

After  the  proper  expansions  are  substituted  into  the  momentum  equation,  the  divergence 
theorem  can  be  used  to  transform  the  Galerkin  equations  into  a  set  of  ordinary  differential 
equations  for  the  components  of  the  velocity  potential.  The  cosine  and  sine  functions 
completely  separate  into  two  uncoupled  sets.  For  the  cosine  elements,  the  Galerkin 
procedure  gives: 
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(2.26) 


where  the  matrices  are  expressed  in  terms  of  integrals  on  the  plane  of  the  rotor. 
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A  similar  set  can  be  written  for  the  sine  components.  Note  that  the  equations  have  been 
partitioned  into  subscripts  “o”  for  terms  with  m+n  odd  and  “e”  for  terms  with  m+n  even. 
For  cosine  terms,  m— 0,1, 2, 3  .  .  .  and  for  sine  terms,  m— 1,2,3...  In  either  case,  n—m,  m-‘r2, 
m+4  orm+1,  m—3,  m+5  depending  on  whether  the  subscript  of  the  partition  is  “d’  or  “o.” 


The  alternative  forms  in  Equations  (2.27-2.28)  (along  with  identities  for  the  Legendre 
functions)  allow  the  matrices  to  be  reduced  to  integrals  on  the  rotor  disk  which  can  be 
expressed  in  closed  form.  Thus,  the  [D]  matrix  can  be  expressed  as  below  for  either  cosine 
or  sine; 
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matrix  can  be  expressed  as: 
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X  =  tan(j/2),  /  =  min(r,/^) 
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It  should  be  noted  that,  because  of  the  logarithmic  term  in  the  m—n—0  element,  the 
integral  for  /"is  formally  infinite  for  the  case  /^  =  r=  »=y  =  0.  However,  following  the  lead 
of  Ref.  [4],  we  express  the  infinity  as  a  finite  series  that  only  approaches  infinity  as  the 
number  of  terms  in  the  expansion  approaches  infinity,  but  that  is  finite  for  the  case  of  a 
truncated  approximation.  This  allows  a  solution  with  a  finite  number  of  terms  that  provides 
formal  convergence  to  the  exact  answer, 

c=f4i;-l4 

\n  „^^n)  2 

where  //  is  the  largest  subscript  appearing  in  the  equations. 

Thus,  the  above  equations  give  a  closed-form  representation  of  the  dynamic  equations 
for  the  velocity  potential  expansions  in  terms  of  prime  potentials.  Therefore,  Equation 
(2.26)  and  its  sine  counterpart  are  a  complete  set  of  equations  for  rotor  induced  flow. 

2.5  Equations  in  Terms  of  Derived  Potentials 

The  next  step  is  to  make  a  change  of  variable  from  a  series  in  terms  of  the  prime 
potentials.  Equation  (2.12),  to  a  series  in  terms  of  derived  potentials.  Equation  (2.18).  The 
change  of  variable  is  effected  through  a  Galerkin  procedure  that  relates  the  two  sets  of  basis 
functions  over  the  upper  hemisphere.  In  particular,  the  total  velocity  potential  is  written  in 
the  two  bases  as  below. 
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{m)  ^  r \T/ mc~\  {  ^  f /T\ fn  ,  /T\ ^  1 

a„]  \^n  \  =[a„]  [^„  ®«+l  +  <;„  } 


(2.41) 


Then,  the  equations  are  dotted  with  the  gradient  of  each  potential  function  and  integrated 
over  the  upper  hemisphere,  as  with  the  momentum  equations.  The  result  is  a  transform 
between  the  two  bases. 

where 


(2.43) 


A  similar  transform  exists  between  b  and  b  that  involves  and  .  Note  that,  in 

axial  flow,  the  transformation  is  identity  since  the  two  potential  functions  and  [^L,J 


are  identical  in  that  case. 


With  the  appropriate  change  of  variable,  the  equation  of  motion  in  terms  of  the  new 
variables  is: 

R[^r]hA+v^[D■][^;J'[M■]{^:]=[D■]{T:■]  (2.44) 

The  version  in  Equation  (2.44)  allows  a  closed-form  expression  for  all  components  of  the 
velocity  field  in  the  entire  upper  hemisphere. 

Interestingly,  when  the  terms  with  m+n  even  are  truncated  from  Equation  (2.44) — and 
when  one  limits  oneself  only  to  the  normal  component  of  flow  at  the  disk — one  recovers 
one  of  the  two  alternated  forms  of  the  He  model  of  Ref.  4,  namely  the  form  with  the 
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velocities  expanded  in  the  P”  [v)  .  The  other  form  of  the  He  equations — the  one  with 

velocity  expanded  in  terms  of  {Xjv^P”  {y^  — is  easily  recovered  by  a  change  of  variable. 

Similarly,  the  Pitt  model  of  Ref.  1  is  recoverable  from  the  second  form  of  the  He  model 
when  it  is  tmncated  to  only  three  states.  Thus,  the  general  form  of  Equation  (2.44)  explicitly 
includes  all  previous  models. 

Although  the  derived  potentials  have  the  advantage  of  being  in  closed  form  everywhere 
in  the  flow  field,  convergence  with  the  derived  potentials  is  slower  at  the  lower  frequencies 
than  is  convergence  with  the  prime  potentials.  Ref.  [4] .  On  the  other  hand,  the  convergence 
of  velocities  with  the  prime  potentials  is  slower  at  the  higher  frequencies  than  is  convergence 
with  derived  potentials.  For  this  cause,  it  seems  reasonable  to  formulate  results  in  terms  of  a 
blending  of  the  two  types  of  potentials,  as  described  below. 

Once  the  equations  are  solved  (in  either  basis).  Equation  (2.42)  shows  how  the  resultant 
states  can  be  transformed  into  either  basis  for  use  in  obtaining  the  flow  field.  We  now 
define  {/}  and  {/}  as  the  blending  vectors  for  the  T  and  T  bases. 

{/l  +  jv'p}'  {/}  (2.45) 

where  the  total  inflow  must  always  add  up  to  the  original  solution, 

{/}  +  [h]"[m^]{/}  =  {4  (2.46) 

[m^]"  [h] {r}  +  {/}  =  [M^y  =  {a}  (2.47) 
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and  where  Equation  (2.42)  has  been  used  for  the  transform.  The  blending  of  the  two  bases 
is  based  on  the  relative  amount  of  unsteady  flow  in  the  dynamic  response.  Thus,  for  low- 
frequency  motion,  it  is  desired  to  have  primarily  the  'F  basis;  whereas  for  high-frequency 
components,  it  is  desired  to  have  more  of  the  T  basis.  It  follows  that  a  logical  manner  to 
accomplish  this  blending  is  to  base  the  relative  proportion  of  T  and  T  on  the  relative 
magnitude  of  the  acceleration  and  velocity  terms  in  the  dynamic  equations. 

To  this  end,  the  acceleration  component  of  flow  field  is  defined  as  follows: 

W=-^[vf'T']W=W-^{r}  (2.48) 

^00  ^  ^00 

Although  the  above  is  only  written  for  the  cosine  component,  one  can  easily  extend  this 
definition  to  include  a  vector  of  both  sine  and  cosine  terms. 

Under  this  definition,  the  magnitude  of  the  blended  vectors  can  be  rationally  defined  as: 

(2.49) 

|(7| 

{/}  (2.50) 

One  can  see  that,  when  the  frequency  goes  to  zero,  acceleration  term  »=0;  and  the  vector 
{/}  =  while  the  vector  |/j  =  0  giving  a  solution  in  terms  of  prime  potentials.  On  the 
other  hand,  at  high  frequencies,  for  which  \a\  is  small,  the  above  formulation  gives  {/}  =  0 
making  the  velocity  expressed  entirely  in  terms  of  which  implies  a  velocity  field  in  terms 
of  derived  potentials.  For  free  vibrations  — in  which  r  =  0 —  {u)  and  [a)  are  equal  such 
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that  the  velocity  field  consequently  is  equally  split  between  prime  potentials  and  derived 
potentials. 


An  implicit  assumption  in  the  above  is  that  the  entire  flow  field  is  at  a  single  frequency. 
In  order  to  accommodate  cases  with  mixed  frequencies  (such  as  a  steady  lift  upon  which  is 
superimposed  a  4/rev  blade-passage  input),  we  can  make  a  simple  modification  to  Equations 
(2.49-2.50).  In  particular,  we  can  treat  the  blending  factors  (on  the  left  of  each  equation)  as 
diagonal  matrices  such  that  each  diagonal  term  (and  each  component  of  flow)  can  have  its 
own  blending  based  on  its  own  component  of  acceleration.  For  example,  each  harmonic 
partition  (i.e.,  superscript,  ni)  could  have  a  factor  based  on  the  Euclidean  length  of  that 
particular  /^-partition: 


/*  = 


{/■}'{< 


a  \  Ka  \  Ku  \ 


(2.51) 


This  allows  the  zero-harmonic  flow  to  be  partitioned  differently  from  (for  example)  the 
four-harmonic  flow. 


2.6  Effect  of  Lift  Tilt 

Makinen  (Ref.  10)  showed  that  the  inflow  model  can  indeed  include  the  effect  of  tilt  on 
the  lift  provided  that  a  correction  is  applied  for  the  swirl  kinetic  energy.  Thus,  the  added 
energy  is  added  to  the  mass  matrix;  and  the  resultant  induced  flow  is  assumed  parallel  to  the 
tilted  lift  vectors. 
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To  be  precise,  for  the  He  model  the  apparent  mass  matrix  (diagonal),  must  be 

replaced  to  include  the  effect  of  the  wake  swirl.  There  are  different  swirl  corrections  that 
can  be  applied,  but  from  Ref.  1 0  the  following  correction  gives  the  best  results  in  axial  flow. 


is  substituted  by: 


where  k  —2.2,  ^  is  the  number  of  blades,  A  is  the  total  inflow,  and  is  the 

transformation  matrix.  It  should  be  noted  that  for  an  actuator  disk  (no  lift  dlt),  k  is  set  to 
zero. 


derived  potentials.  For  the  derived  potentials,  the  modified  mass  matrix  is  found  by 
replacing  the  upper  left  partition  of  the  matrix  [D]  by  the  inverse  of  Eq.  (2.54). 
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2.7  Velocity  in  Lower  Hemisphere 


The  development  on  previous  sections  is  a  Galerkin  method  that  gives  convergence  on 
velocity  in  the  upper  hemisphere  only.  The  next  issue  is  how  to  extend  that  convergence  to 
the  lower  hemisphere.  To  begin  that  development,  it  is  useful  to  consider  the  case  of  axial 
flow,  for  which  closed-form  solutions  can  be  found  for  the  ^-component  of  velocity.  For 
purposes  here,  it  is  sufficient  to  write  the  closed-form  solution  to  Equation  (2.6)  in  axial  flow 
for  a  step  response  in  pressure  on  either  side  of  the  disk.  The  pressure  consequently 
becomes: 

PiZ,0  =  PizHt)  (2.58) 

The  resultant,  closed-form  solution  for  the  ^-component  of  velocity  is  given  by: 
outside  wake:  pV^p  (^,  /)  =  P(^  -  Vj)  -  P(^) 

(2.59) 

inside  wake:  (^,  /)  =  2P(0^ )  -t  P(^  -  Yj)  -  P(^) 

where  the  region  inside  of  the  wake  is  defined  by  T  <  1  and  0  <  ^  <  .  Note  that  the 

above  solution  satisfies  the  momentum  equation,  the  initial  conditions,  and  the  velocity  jump 
across  the  disk. 

pV„h(0")-^,(0-)]  =  P(0")  +  P(0-)  (2.60) 

The  first  term  P(^  is  the  particular  solution  to  the  applied  pressure  field,  and  the  second 
term  Pf^-FG^  is  the  homogeneous  solution  that  satisfies  the  initial  conditions.  One  can 
think  of  the  first  term  as  the  velocity  due  to  the  pressure  drop  at  the  disk  and  the  second 
term  as  the  velocity  due  to  the  shed  vorticity  (like  a  start-up  vortex)  that  is  being  convected 
downstream.  This  solution  satisfies  the  differential  equation  and  initial  conditions,  but  it 
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does  not  satisfy  the  continuity  condition  across  the  disk  (or  across  the  start-up  region). 
Thus,  in  the  wake,  the  extra  term  2P(0^)  — the  pressure  on  the  downstream  side  of  the 
disk — must  be  added  to  enforce  this  condition.  Thus,  the  flow  inside  of  the  wake  is  no 
longer  a  potential  flow.  Interestingly,  for  pressure  drops  across  the  disk,  since  Pf'^  is 
discontinuous  across  the  disk,  this  extra  term  must  be  added  to  enforce  continuity.  For  mass 
sources  at  the  disk,  for  which  pressure  is  continuous,  the  identical  term  must  be  added  to 
enforce  the  velocity  drop  due  to  the  extra  flow.  Equation  (2.60).  Thus,  the  same  added  term 
holds  for  aU  pressure  inputs  — lift  or  mass  sources. 


For  example,  for  a  step  input  in  pressure  drop  (and  after  a  sufficient  amount  of  time), 
comparison  of  the  “outside  wake”  and  “inside  wake”  solutions  shows  that  the  velocity  is 
continuous  cross  the  disk  and  is  given  from  the  relationship  pY^v  =  — P(0  )  =  +P(0^) .  This 
velocity  then  expands  downstream  to  2P(0^)  as  the  pressure  drops,  a  well-known  property  of 
rotor  wakes  due  to  the  Bernoulli  equation.  Further  downstream,  as  one  approaches  the 
start-up  system  (^  =  the  velocity  diminishes  back  down  to  P(0^)',  and  then  (outside  of 
the  wake)  it  further  diminishes  down  to  zero  in  the  far  field.  It  is  the  added  term  in  the 
“inside  wake”  solution  that  allows  for  this  correct  description.  This  concept  is  easily 
extended  to  the  general  case  through  the  use  of  a  convolution  integral  involving  the  step 
response.  Thus,  for  a  general  pressure  input  P(%^,  the  difference  between  the  “outside 
wake”  and  “inside  wake”  ^-velocity  fields  is: 


v..=  p\ 

^  inside 


^  ^  \  outside 


+  2P(0",/-^/VJ 


(2.61) 
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One  can  see  that  the  added  term  is  proportional  to  the  downstream  pressure  that  was  on  the 
disk  at  a  past  time  (delayed  by  how  long  it  took  the  wake  to  travel  from  the  rotor  to  that 
particular  ^location). 

From  the  above  development,  one  can  formulate  an  approximation  to  the  velocity  field 
below  the  plane  of  the  rotor.  For  points  outside  of  the  wake.  Equation  (2.59)  shows  that  the 
velocity  is  basically  the  analytic  continuation  of  the  velocity  field  above  the  disk.  Thus,  one 
can  simply  continue  the  potential  functions  into  the  lower  hemisphere  and  find  that  velocity 
field  there.  Admittedly,  the  velocity  expansion  in  the  lower  plane  will  converge  more  slowly 
than  in  the  upper  hemisphere  with  the  convergence  deteriorating  as  one  approaches  the 
downstream  end  of  the  wake.  Flowever,  the  velocity  is  decaying  in  that  region  so  that  errors 
are  minimized.  For  velocities  within  the  wake,  one  can  take  the  extended  lower  hemisphere 
solution  and  add  the  term  in  Equation  (2.61)  in  order  to  obtain  an  approximation  for  within 
the  wake.  One  again,  it  is  emphasized  that  this  is  only  an  approximation  in  that  it  neglects 
the  shed  “end  wakes”  for  each  oscillation.  For  example,  for  a  time-varying,  square-wave 
pressure,  the  above  procedure  will  give  the  correct  square-wave  behavior  in  the  wake  only  at 
low  frequency.  Flowever,  as  frequency  increases,  and  as  the  shape  of  the  velocity  wave 
begins  to  warp  with  respect  to  a  square  wave,  the  above  procedure  will  be  slow  to  converge 
on  that  warping. 

To  put  the  above  in  more  concrete  terms,  for  a  rotor  with  lift  concentrated  at  a  finite 
number  of  concentrated  blades,  the  inflow  between  the  vortex  sheets  should  be  modeled 
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fairly  accurately  by  the  above;  but  the  discontinuity  at  the  vortex  sheets  will  be  slow  to 
converge.  Although  this  is  not  a  perfect  solution,  it  seems  to  be  a  practical  solution. 

Furthermore,  the  above  approach  is  probably  the  best  that  can  be  done  with  closed-form 
velocity  potentials  and  a  Galerkin  method.  For  skewed  flow,  the  velocity  below  the  plane  of 
the  disk  is  treated  in  entirely  the  same  way  with  the  accompanying  approximations.  Of 
course,  in  the  limit  as  wake  skew  angle  approaches  90°,  the  methodology  breaks  down  as  the 
wake  goes  to  2ero  thickness  and  does  not  die  out  (even  in  the  upper  hemisphere).  Flowever, 
it  should  be  pointed  out  that,  even  for  that  case,  the  flow  on  the  disk  stiU  converges.  It  is 
only  off-disk  flow  that  cannot  converge  for  90°  skew  angle.  Refs.  [4-8].  It  should  be  noted, 
however,  that  this  approach  does  give  the  exact  result  for  the  flow  at  all  wake  skew  angles 
(even  90°)  for  the  case  of  steady  flow  in  the  wake.  Thus,  it  is  perfectly  suited  to  issues  of 
wake  interactions  with  lifting  surfaces. 

In  essence,  the  above  development  shows  that  one  would  require  more  states  to  model 
the  flow  in  the  wake  than  one  would  need  to  model  the  flow  outside  of  the  wake.  This 
conclusion  follows  from  the  fact  that,  in  a  state-space  approach,  a  pure  time  delay  must  be 
treated  by  extra  states;  and  the  above  development  involves  a  time  delay  in  the  wake. 
Therefore,  for  future  work  in  this  area,  attention  would  need  to  be  given  to  the 
incorporation  of  additional  vortex  states  within  the  wake. 
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2.8  Nonlinearities 


The  above  development  is  for  flows  that  are  linearized  about  Although  this 

assumption  is  quite  good  in  cruise  conditions,  it  breaks  down  in  the  limit  of  a  helicopter  in 
hover.  At  the  hover  condition,  goes  to  0;  and  the  equations  become  singular.  From  the 
very  inception  of  dynamic  wake  models,  this  was  recognized.  Fortunately,  there  is  a  fairly 
simple  way  to  make  the  equations  nonlinear  in  such  a  way  that  they  agree  exacdy  with 
momentum  theory  for  the  hover  condition  and  then  transition  to  the  linearized  case  of  Ref. 
[1].  To  begin  this  development,  it  is  useful  to  replace  the  term  with  the  total  average 
flow  at  the  rotor. 

K,.  =  sin"  j  +  iv^  cos X  +  \)  (2-62) 

Flowever,  this  definition  is  not  sufficient  to  create  a  model  that  will  exactly  linearize  to 
the  correct  equation.  For  the  correct  linearization,  one  must  consider  the  perturbation  of 
the  complete  nonlinear  equations  which  involves  perturbation  of  the  parameter  defined  as 
the  mass-flow  parameter,  V. 

d  /  xKfsin^  J  +  (KoCosj  +  fJ(14cosj  +  2fJ 
dv}'^  '‘d  IZ,.  ^  ^ 

The  mass  flow  parameter  allows  the  creation  of  a  set  of  nonlinear  equations  that  will  agree 
with  momentum  theory  in  hover  and  also  linearize  to  the  correct  equations  derived  above. 
In  particular,  one  defines  a  diagonal,  {V\  matrix. 
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The  above  equations  are  nonlinear  about  the  elliptical  pressure  distribution  r”  and  thus 
must  use  only  the  average  inflow  due  to  that  lift  in  determining  the  U  to  go  into  the  Vj  and 
V  parameters. 

[LY[M]{a]  (2.61) 

The  result  is  that  the  equations  linearize  to  the  equivalent  of  the  original  linear  equations 
but  with  l/oo  of  those  equations  replaced  by  V  at  all  instances.  This  then  reduces  to  the  Pitt 
model  of  Ref.  [1]  when  only  three  states  are  taken. 

As  a  bi-product,  one  can  see  that  (due  to  continuity)  the  nonlinear  version  implies  a 
contracting  of  the  rotor  wake  downstream  as  the  induced  velocity  increases.  This  can  readily 
be  handled  by  redefining  the  “wake”  region  has  having  a  smaller  radius  than  the  rotor.  In 


particular,  if  one  computes  the  average  velocity  at  any  circular  cross-section  of  the  wake  (the 
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cut  made  parallel  to  the  ^-plane)  and  then  calls  that  average  velocity  t'  (^)  ,  then  we  can  write 
the  equivalent  radius  as: 


vRy 


m+K 

m+K 


(2.62) 


where  ^  will  be  negative  in  the  wake.  This  will  ensure  continuity.  The  flow  outside  of  this 
new  wake  deflnition  will  remain  the  out-of-wake  solution 


The  above  is  a  nonlinear  model  for  velocity  that  follows  Refs.  [1-4]  in  how  the 
nonknearities  due  to  hover  are  treated  in  the  velocity  field.  However,  none  of  those 
references  discusses  how  the  nonlinear  pressure  fleld  differs  from  the  linearized  one.  The 
difference  is  clearly  seen  for  the  lifting  case  in  hover  for  which  the  linear  theory  predicts 
equal  and  opposite  pressures  on  either  side  of  the  disc  — giving  a  total  pressure  increase  of 

2P^0^^  =  — 2P^0  Momentum  theory,  on  the  other  hand,  shows  that  the  positive 

pressure  below  the  disk  is  y4  of  the  total  pressure  rise  while  the  negative  pressure  above  the 
disk  is  y4  of  the  total.  The  reconciliation  of  this  anomaly  is  easily  found  by  substitution  of 
the  nonlinear  velocity  fleld  into  Bernoulli’s  equation.  With  this  procedure,  the  nonlinear 
pressure  in  the  flow  field  can  be  written  in  terms  of  the  linearized  pressures  (written  in  terms 
of  z)  as  follows. 


Outside  of  wake: 


P 

P 


P  1  ^  ^ 

—  --v-v-v^V^ 

P  ^ 


(2.69) 


Inside  of  wake: 


P 

P 


P,  2P(0^)  1  -  ^ 

—  + - 

Po  P  2 


(2.70) 
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where  Pq  is  the  total  atmospheric  pressure  in  the  far  field,  is  the  free-stream  velocity, 
is  the  pressure  (above  atmospheric)  applied  below  the  disk  (whether  from  lift 
differential  or  mass  sources  at  rotor),  p  is  the  nonlinear  induced  velocity  from  the  nonlinear 
inflow  theory  described  above,  and  is  the  component  of  p  along  the  negative  ^  direction 


(the  free-stream  direction  of  For  the  perturbation  (i.e.,  linearized)  model,  if  follows 

that: 


Outside  of  wake: 


P 


(2.71) 


Inside  of  wake: 


P 

P 


2P{0^) 

P 


-  p-p-p^V^ 


(2.72) 


where  P  is  the  perturbation  pressure,  p  is  the  perturbation  induced  velocity,  and  p^  is  the 
perturbation  velocity  in  the  free-stream  direction,  and  ¥  is  the  steady  induced  flow  about 
which  the  equations  are  being  perturbed.  This,  then,  completes  a  non-linear  model  and  a 
perturbation  nonlinear  model  of  the  inflow  theory. 


2.9  Results  from  previous  models 


Many  results  from  this  unified  approach  are  available  (Refs.  1-14)  including  comparisons 
with  experimental  data  and  with  other  theories.  The  results  present  velocity  components  in 
the  plane  of  the  rotor  disk  for  the  linear  equations  with  simple  harmonic  pressure  or  mass- 
source  inputs.  In  all  results,  ten  odd  functions  and  ten  even  functions  are  used  in  the 
Galerkin  computations.  For  the  case  of  frequency  response,  there  is  an  exact  numerical 
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solution  based  on  a  convolution  integral  along  streamlines  from  far  upstream  to  the  point  at 
which  induced  flow  is  desired.  Thus,  the  comparisons  to  follow  are  validation  that  the 
potential  functions  and  the  Galerkin-based  equations  are  correct. 

A  complete  dynamic  wake  model  has  been  derived  from  a  Galerkin  procedure  applied  to 
the  potential  flow  equations.  The  resultant  velocity  potentials  give  converged  solutions  both 
above  the  disk,  on  the  disk,  below  the  disk  and  in  the  wake.  Previous  inflow  models  such  as 
momentum  theory,  Pitt  model,  and  the  He  model  are  contained  as  special  cases. 
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Chapter  3 


Development  of  the  Solution  Approach 


In  an  attempt  to  understand  how  the  tilt,  the  number  of  blades  and  the  geometry  of  the 
blades  in  a  real  rotor  affect  the  induced  power  the  present  research  will  expand  He’s 
generalized  dynamic  wake  model  to  study  minimum  induced  power  for  a  rotor  under  a 
variety  of  conditions.  The  scope  of  such  a  study  is  so  broad  that  the  use  of  large, 
comprehensive  codes  is  prohibitive  for  these  purposes.  On  the  other  hand,  finite-state  wake 
models  are  ideally  suited  for  such  task.  These  models  expand  both  the  pressure  field  and  the 
velocity  field  in  orthogonal  expansion  functions.  Therefore,  the  computation  of  induced 
power  (the  dot  product  of  thrust  and  induced  flow)  simplifies  nicely  into  a  quadratic  cost 
function  that  allows  classical  optimization  to  be  used  for  the  minimum  power  under  a  variety 
of  constraints.  Thus,  it  is  anticipated  that  such  an  approach  can  yield  insight  into  this  issue. 
As  a  preliminary  step  in  such  an  endeavor,  this  present  dissertation  looks  at  the  induced 
power  of  a  non-ideal  lifting  rotor  in  axial  flow  to  verify  that  dynamic  wake  models  can 
indeed  compute  the  proper  induced  power. 

Since  theory  and  experiment  agree  with  simple  momentum  approaches  for  power  in  axial 
flow,  such  conditions  provide  the  ideal  test  bed  to  verify  that  this  optimization  approach  is 
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viable.  The  remaining  of  the  research  efforts  will  then  concentrate  on  results  for  induced 
power  in  forward  flight,  using  the  formulation  already  derived  and  included  on  this 
dissertation. 


3.1  Induced  Power  Derivation  from  He-Peters 


He  [4]  developed  an  unsteady  induced-flow  theory  to  be  used  in  stability,  vibration, 
control,  and  aeroelastic  studies.  The  theory  is  based  on  an  acceleration  potential  for  an 
actuator  disk.  The  induced  flow,  w,  is  expressed  in  a  polynomial  distribution  (proportional 
to  Legendre  functions)  radially  and  in  terms  of  a  Fourier  series  azimuthally.  The  way  the 
induced  flow  is  set  up  allows  for  all  harmonics  and  describes  the  induced  flow  for  any  radial 
position. 

The  theory  provides  the  pressures  on  the  rotor  disk  as  a  Fourier  expansion.  As  more 
harmonics  are  added,  that  pressure  converges  to  the  lift  concentrations  on  the  blade  and  to 
zero  lift  off  the  blade.  One  of  the  issues  to  be  addressed  in  this  present  work  is  whether  or 
not  such  an  approach  can  give  adequate  convergence  to  induced  power.  The  form  of  this 

pressure  expansion  is  as  follows: 

00  00 

m=^  n=m-\-\,m-V'h,... 

The  pressure  at  the  rotor  disk  is  obtained  by  the  difference  between  the  pressure  above  and 
the  pressure  below  the  disk. 
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ou  ou 

AP{r,i^,0  =  Z  P:(y)[C(T)cos(mi/^)  +  C(T)sin(mi/^)]  (3.2) 


m=0  fi=m+l,m+3,... 


The  power,  can  be  expressed  as: 


^  =  II AP  (^ip  +  vjxdx^ii// 


The  Peters-He  model  also  sets  out  the  velocity  field  normal  to  the  rotor  disk  in  terms  of  the 
same  Legendre  Functions  and  variables,  as  given  below. 


00  00 

'j(r,i//J)  =  QR'^  Z  €ir)[_a”{J)(:o?,{my/)  + p”{J)sin{my/)\  (3.4) 


=0  «=/^+l,«?+3,... 


€(r)  =  ^P:(y)  (3.5) 

where  /  =Q/,  and  ^”(T)are  a  complete  set  of  functions  that  arise  from  the  solution  to 
Laplace’s  equation  in  ellipsoidal  coordinates. 

The  form  of  the  functions  is: 

=  +  Z  ^  - 7TT7  (3-6) 

?=».»+2....  (^  - 1^)  !!(^  + 1^) !!(«  -  ^  - 1) !! 


where  H"  = 


substitution  of  the  induced  flow  and  the  pressure  at  the  disk  yields: 


5  =  £  A  +  '^-Pj(v)(a'^cos(ri//)  +  /?-sin(ri//)J 


P”  (v')(r*  cos{m\i/)  +  r”  sin^my/)^  vdvdy/ 
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(3.8) 


5  = 


j  j 


+7rpQ.^R^  m{[  ^ j  ^  j  ^  n  j  ^  n  ^ 


The  functions  P”  are  Legendre  Functions  of  the  first  kind,  and  V  is  related  to  the  radial 


position 


be  r  =  "s/l^ 


V  . 


The  equations  that  relate  the  pressure  coefficients  in  the  pressure  expansion 


to 


the  velocity  coefficients  are  derived  from  the  momentum  equation  of  potential 


flow. 


[k;  ]  {a;} + k  [£'  ]"  {<}  =  he } 


(3.9) 


[V']iA'f+'^[r'] ‘(A'}  =  hc'} 


(3.10) 


where  (  )  =  — ,  IT  =  ^  ^  — ,  A  is  the  total  inflow,  u  is  the  advance  ratio,  IT  is  the 


2 

flow  parameter,  and  K”  is  diagonal;  K”  =  —  H"  .  The  L  cosine  and  sine 

n  ^ 


matrices  are 


given  in  closed  form  in  terms  of  the  wake  skew  angle,  . 


Jn 


n+ j—2r 

(-1)~  2^{2n  +  l){2j+l) 

^H;h;  (7  +  «)(7  +  «  +  2)[(7-«)"-1] 


for  r  +  m  even 


jn 


n 


sign{r  —  m) 


2^H:h;  j(2»  +  \X2j+\) 


for  r  +  m  odd,y— »±  / 
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_  Q 

jn 


for  r  -‘r  m  odd,y^»±  / 


(3.11) 


^jn  J  “  ^  jn 

[LjJ  =rxi*'''i  -(-!)' 


(3.12) 


where  /—  min(r,m),  X  —  tan|)(/2|.  The  forcing  functions,  r” ,  are  given  in  terms  of  the 
blade  loading,  L^.  These  are  related  to  the  pressure  and  the  optimization  of  the  induced 
power  wiU  yield  an  optimal  pressure  distribution  for  a  given  thrust  coefficient. 


C=-S  Jo 

^  q=\  /X2  K 


^  0.  it 


(3.13) 


For  the  Peters-He  model  in  its  actuator-disk  form,  we  have  a  skewed  wake  as  shown  in 


Figure  3.1  below. 
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Inflow,  w 


Figure  3.1:  Normali2ed  velocity  components. 

The  power  ^  can  be  found  from  Equation  (3.3)  to  be: 

s=|J  xdxdy/ 

A 

Notice  that  the  power  does  not  depend  on  the  velocity  component  u.  The  pressure  at  the 
disk  (as  shown  in  Equation  (3.2))  is: 


AP{r,\i/,t)  =  pQ.^R^Yj  X  ^/(^)[C(Ocos(/^^^)  +  r"(/)sin(^^^)] 

w=0  n=m+\,m+Z,... 


The  advance  ratio  is  a  non-dimensional  parameter  used  for  a  rotor  in  forward  flight  that 
is  the  ratio  of  the  axial  component  of  the  incoming  velocity  and  the  rotational  velocity  of  the 
rotor,  shown  in  Equation  (3.14)  and  sketched  in  Figure  3.1.  Some  other  significant 
parameters  are  defined  for  convenience. 

.a  •  ^ 

Advance  ratio  U  = - 

QR 

Climb  rate  A.  = - 

QR 


(3.15) 


X 

Non-dimensional  radial  position  ^ 

(3.16) 

—  w 

Non-dimensional  inflow  K= - 

QR 

(3.17) 

Note  that  rdr  — vdv,  and  X  and  [t  are  constant. 


The  inflow  is  given  by  Equation  (3.4).  Introducing  the  definitions  for  pressure  change,  p 
from  the  cMmb  rate,  and  the  induced  velocity  we  obtain  the  expression  for  the  power  where 
the  cUmb  rate.  A,  is  a  constant.  The  normalized  Legendre  function  is  =  Vsv  by  definition. 
Introduction  of  it  into  the  expression  for  the  power  yields  Equation  (3.18).  By  the  use  of 
the  simple  relationship  between  power  and  power  coefficient,  we  obtain  the  following. 


n 


1 


;L+y-p 


j  ^  j 


\^y:ydv 


m,n  j 


-P”’a” 

J  J 

V 


Pm  ^mc  . 

T  + 

n  n 


1  - 

^  J  '  j 
V 


(3.18) 


pyr  \ydv 


Q  = 


5 


■  =  2 


S 


+«, 


Oc 


«=3,5,.. 


aV:  + 


I  z 

.w=l,2,3,...  J=m+\,m+'i,... 


^  m  _mc 

a  , 

J  J 


^an,  m 
J  J 


■'  (3.19) 


and 


~ — Jo  «’Zt,'(t')[c  co^i^mxj/')  +  r"  sin(;w^)J  vdvdy/  (3.20) 


A  simple  check  using  m  =  0  only,  and  n—\  only,  provides  the  common  expression  that 
shows  the  Peters-He  model  agrees  with  the  induced  power  from  Momentum  Theory.  The 
lift  and  pressure  coefficients  for  this  special  case  are  shown  in  Equations  (3.21)  through 
(3.23). 
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(3.21) 


;'=^c 

2  ' 


(3.22) 


Cp  -  +  -C-p{A  +  v) 


(3.23) 


Equations  (3.19)  and  (3.20)  provide  the  framework  for  a  classical,  quadratic  optimization  of 
power  subject  to  constant  thmst. 


3.2  Induced  Power  Theorem 


A  Rotor  Induced-Power  Theorem  is  used  to  verify  the  approach  of  this  work.  Let  a  rotor, 
Figure  3.2,  be  moving  along  an  arbitrary,  straight  path  through  still  air  with  a  velocity  W.  Let 
X  be  the  angle  between  the  flight  path  and  a  vertical  to  the  rotor. 


Figure  3.2:  Inflow  velocity  components  for  a  moving  rotor. 
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It  follows  that  the  inplane  component  of  air  velocity  as  seen  by  the  rotor  is  U  —  IFsinj 
and  the  normal  component  is  V  =  Ihcosj. 

Let  the  rotor  loading  perpendicular  to  rotor  plane  be  called  T  and  the  rotor  load  in  plane 
be  called  H,  each  with  a  positive  sense  in  the  direction  of  the  flight  path  (i.e.,  opposite  to  V 
and  U).  Let  the  blades  in  the  rotor  disk  be  rotating  counter-clockwise  at  angular  velocity  Q, 
when  looking  down  on  the  rotor,  and  let  ^  be  the  azimuth  angle  of  a  blade  as  measured  from 
aft,  tj)  =  Qt.  Let  a  generic  point  on  the  blade  be  a  radial  distance  x  from  the  center  of 
rotation  as  shown  on  Figure  3.3. 


Figure  3.3:  Rotor  blade  generic  position  and  rotating  angle. 


Let  (p  be  the  inflow  angle  as  seen  in  the  local  blade  system.  Figures  3.3  and  3.4.  In  that 
system,  let  dh  be  the  incremental  local  lift  per  unit  length  (perpendicular  to  the  total  inflow). 
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let  dD  be  the  incremental  local  induced  drag  per  unit  length,  and  let  cfT  be  the  incremental 
thmst.  Let  w  be  the  induced  flow,  opposite  to  -L,  Figure  3.4. 

dT 


Figure  3.4:  Geometry  of  the  forces  on  the  blade. 


dT  =  dT  CO?,  (p 

(3.24) 

dD  =  dT  sin  (p 

(3.25) 

dH  =  —dD  sin  if/ 

(3.26) 

Figure  3.5,  taken  from  the  work  of  Glauert  [15],  shows  the  geometry  of  the  flow  in  the 
blade  coordinate  system.  The  relative  flow  due  to  rotor  motion  alone  is  Qx  +Usm(y/)  in  the 
rotor  plane  and  V  perpendicular  to  that  plane.  The  induced  flow  n>  must  be  parallel  to  the 
lift,  so  it  is  added  vectoraUy  at  the  angle  tp  as  shown.  The  resultant  total  inflow  (due  to  rotor 
motion  and  due  to  induced  flow)  must  be  perpendicular  to  the  local  incremental  Uft,  due  to 
circulation  considerations.  Therefore,  w  can  be  considered  perpendicular  to  the  total  flow 
vector.  The  resultant  relationship  gives  rise  to  the  geometry  in  the  flgure  and  to  the 
following  identities: 
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w 


WCO'iCp-^V 

tan  (p  = - ^ ^ - 

Qx  —  w  sin  (p  +  U  sin  y/ 


V  + 

cos  (p 

Qx  +  U  sin  Xj/ 


(3.27) 


V  cos(p  +  w 

sin^  =  tan^cos^  = - ^ -  (3.28) 

Qx  +  U  sin  Xj/ 

The  above  can  be  used  to  transform  the  induced  power  equations. 

=p^-TV-HU  (2.1) 

Pj  =  LQx  sin  ^  —  _L  cos  +  _L  sin  ^  sin  Xj/U  (3.29) 

Pj  =  _L(Qx  +  U  sinxj/)  sin  (p  —  _L1T  cos  (p  (3.30) 

(where  the  differentials  are  omitted  for  clarity). 


.  V  cos  (p  +  w 

But  sin  (p  = - ^ - .  Therefore 

Qv  +  U  sin  Xj/ 


Pj  =  LV  cos  (p  + 1^(0  —  l^V  cos  (p  =  Yjv 


(3.31) 
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The  induced  power  is,  then: 


Pj  =  Lip 


(work  done  by  L  on  2p)  (3.32) 


Thus,  the  incremental  induced  power  can  be  found  from  the  integral  of  the  dot  product 
of  the  local  lift  and  local  induced  flow,  which  is  the  work  done  on  the  flow  field.  The  above 
theorem  is,  strictly-speaking,  exacdy  tme  only  for  axial  flow  because  of  the  assumption  that 
local  lift  is  parallel  to  local  induced  flow.  On  the  other  hand,  that  assumption  is  less  and  less 
important  as  one  transitions  away  from  hover.  Furthermore,  it  is  exactly  true  that  the  work 
done  on  the  flow  field  will  equal  the  induced  power.  Therefore,  Equation  (3.32)  seems  a 
valid  approach  to  computing  the  induced  power  from  a  dynamic  wake  model. 


3.3  General  Optimization 


A  general  classical  quadratic  optimization  problem  is  stated  as  follows: 

Minimize  {x}*  [M]{x}  subject  to  {x}  =  ^  (given).  Use  of  Lagrange’s  multiplier  to 
include  the  constraint  leads  to  the  cost  function: 

J  =  (3-33) 

where  A  is  the  Lagrange  multiplier.  Optimizing,  we  obtain  that  for  the  change  of  the 
functional  to  be  zero 


SJ  =  {^x} 


=  0 


(3.34) 
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X  = 


f}  A 


(3.35) 


Notice  that  the  matrix  to  be  inverted  is  the  symmetric  part  of 


The  Lagrange  multiplier  must  be  chosen  such  that: 


A  = 


2W'TM]+Mr 


(3.36) 


The  above  derivation  provides  a  general  guide  for  classical  quadratic  optimizations,  such 
as  the  one  that  is  performed  for  the  present  research.  The  Lagrange  multiplier  will  be 
chosen  as  specified  by  equation  (3.36)  for  our  special  case  of  an  optimum  rotor. 


3.4  General  Figure  of  Merit  Definition 


Figure  of  merit  is  a  measure  of  rotor  efficiency.  It  is  defined  as  the  ratio  of  minimum  or 

ideal  power  required  to  hover  to  the  actual  power  required.  It  compares  actual  rotor 

performance  with  the  performance  of  an  ideal  rotor. 

^  ^  minimum  possible  power  required  to  hover 
actual  power  required  to  hover 
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Note  that  the  ideal  figure  of  merit,  or  the  figure  of  merit  for  an  ideal  rotor,  should  equal 
unity.  The  closer  the  figure  of  merit  is  to  unity  for  a  given  rotor,  the  more  efficient  the  rotor 
is. 

In  this  thesis,  we  extend  this  definition  of  the  figure  of  merit  to  climb  and  to  forward 
flight  by  the  definition 

^  ^  minimum  induced  power  at  any  flight  condition 
actual  induced  power  at  that  flight  condition 
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Chapter  4 


Results 


The  classical  results  for  the  cases  of  axial  flow,  plus  the  case  of  an  actuator  disk  with  an 
infinite  number  of  blades  in  skewed  flow  are  given  in  this  section.  The  finite-state  method 
used  in  this  research  needs  to  be  verified.  For  the  purpose  of  verification,  solutions  are 
obtained  from  classical  assumptions  and  procedures  to  be  compared  to  the  results  obtained 
by  the  current  approach. 


4.1  Classical  Theories:  Axial  and  Skewed  Flows 


4.1.1  Momentum  Theory 


Momentum  theory  utilizes  the  basic  fluid  mechanics  conservation  laws  (Ref.  25)  to 
predict  the  behavior  of  a  rotor,  modeled  as  an  actuator  disk.  Conservation  of  mass, 
momentum  and  energy  equations  yield  results  for  optimum  induced  velocity  and  pressure 
distributions  for  an  actuator  disk  in  hover  or  climbing  and  a  given  thrust  [16]. 
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Consider  an  actuator  disk  of  area  A,  with  a  thrust,  T,  moving  vertically  with  a  velocity  V 
as  shown  in  Figure  4.1.  Consider  a  well  defined  slipstream  that  goes  from  an  area  of  Aj,  far 
upstream,  to  an  area  A2  far  downstream  of  the  actuator  disk  Across  the  disk,  the  velocity 
increases  to  (V  +  v)  due  to  the  change  in  pressure  at  the  location  of  the  disk.  Downstream, 
the  velocity  is  the  addition  of  the  incoming  and  induced  velocities,  (A  +  n).  Applying 
conservation  of  momentum,  we  obtain  that  the  power  is: 

P  =  ^^p{V  +  v)dA  (4.1) 

where  Ap  is  the  pressure  difference  across  the  rotor. 


Applying  Bernoulli’s  equation  upstream  and  downstream  of  the  actuator  disk  along 
streamlines,  the  expression  for  the  pressure  difference  across  the  rotor  is: 


(4.2) 


where  g  is  the  density. 

Therefore  the  power,  P  becomes: 


P  = 


){V  +  v)dA^ 


(4.3) 


Applying  conservation  equations,  we  obtain  the  thmst  to  be: 

T  =  J  p{V  -^-w^wdA^ 


(4.4) 
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The  objective  is  to  obtain  the  minimum  power  distribution  (equation  (4.3))  for  a  given 
thmst  (Equation  (4.4)).  Basically,  this  is  an  optimization  problem  for  a  functional,  P,  subject 
to  a  constraint,  T.  We  would  like  to  obtain  the  induced  velocity  distribution  tp{Y)  that 
minimizes  the  power.  Euler’s  equation  below  can  be  used  to  find  the  distribution  tp{£). 


A 

dr 


V  +  t')  +  /l(yC>(K  +  w') 


7  ^  f  '\ 

—  — -  —  pilVw +  v^-\-  X{^p{V  + 
dt  dr  \  2 


where  A  here  is  the  Lagrange  multiplier. 


Ai 


V+w 

Figure  4.1:  Flow  model  for  an  actuator  disk  in  axial  flow. 

Since  the  expressions  are  independent  of  r,  it  follows  that  the  solution  for  the  minimization 
problem  is  that  w  be  constant.  The  optimum  induced  velocity  distribution  is  therefore  a 
constant  or  uniform  distribution.  If  this  value  is  introduced  into  Equation  (4.2),  the 
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resulting  pressure  distribution  that  produces  minimum  induced  power  is  found  to  be  a 
constant  or  uniform  distribution  also. 


From  momentum  theory  for  uniform  induced  flow  the  thmst  and  power  coefficients  are 
defined  as: 


C,.  =2{ri  +  v)v 
Cp  =2{r}  +  vf  V  =  {r}  +  v)C.p 


(4.6) 

(4.7) 


In  Equations  (4.6)  and  (4.7),  we  normalize  the  velocities  and  power  on  as  follows 

7]  -  V2Cp 


—  V'  _  '/  7^ 

V  =  77  =  ,  Cp  =  — yyf 

C-j- 


.1^'  V772 


(4.8) 


Since  Equations  (4.6)  and  (4.7)  are  the  equations  for  uniform  (which  is  minimum  power),  we 
can  solve  for  the  minimum  possible  Cp  (total  power)  and  Cp—rj  (induced  power) 


Cp=^+  1  + 

2  V 


f  rj'^^ 
\2j 


Cpi  Cp  TJ  ^  +  Jl  + 


f ;;  Y 


(4.9) 


(4.10) 


It  is  helpful  to  rewrite  Equation  (4.10)  by  multiplying  numerator  and  denominator 


?7  (  ^  ■ 

by —  + ,  1  +  —  .  This  gives  an  alternate  form  of  the  minimum  induced  power 

2  V  [2j 


C  = 

PI 


^  L  1 

6/7^ 

2 

— + ji+ 

2  V  I 

(4.11) 
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For  axial  flow  {rj  =0),  the  minimum  Cpj  is  1.0.  This  is  why  the  flgure  of  merit  in  hover 

can  be  defined  as  [Cp  .  Note  that  in  cUmb,  the  minimum  induced  power  is  less  than 

unity  and  approaches  l/?7  for  rj  » 1 .  Thus,  the  figure  of  merit  in  axial  flow  is 

/^Factual  when  Cpj  is  from  Equation  (4.11). 

Comparison  of  momentum  theory  results  with  finite  state  methods  can  be  performed  by 
the  use  of  optimum  pressure  and  induced  velocity  profiles.  Therefore,  the  plot  for  this 
minimum  induced  power  will  not  be  shown.  However,  it  is  important  to  see  the  expression 
for  it  as  it  provides  a  baseline  against  which  to  measure  efficiency,  and  we  can  obtain 
valuable  information  such  as  the  tendency  of  the  efficiency  with  respect  to  the  climb,  the 
shape  of  the  curve,  etc.. 

A  similar  approach  is  followed  to  obtain  the  optimum  induced  velocity  profile  and  the 
optimum  pressure  distribution  for  skewed  flow.  As  in  hover,  the  induced  power  or  power 
loss  can  be  analy2ed  from  the  induced  velocity.  The  difference  is  that,  while  in  forward 
flight,  the  edgewise  velocity  component  of  the  incoming  flow  has  to  be  accounted  for  in  the 
calculations.  The  induced  velocity  is  assumed  to  be  a  constant  in  hover. 

Momentum  theory  shows  that  the  optimum  pressure  distribution  for  an  actuator  disk  with 
an  infinite  number  of  blades  in  skewed  flow  is  a  constant,  whereas  the  inflow  velocity 
distribution  is  no  longer  a  constant. 
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4.1.2  PrandtPs  Approximation 


In  1919,  Betz  and  Prandtl  [18]  combined  momentum  and  blade-element  theory  to  obtain 
a  correction  factor,  k  that  must  be  added  to  the  absolute  performance  calculations  for  a 
rotor  to  account  for  tip  losses.  The  tip  losses  occur  due  to  the  reduction  in  the  thrust  near 
the  blade  tips  when  a  system  with  a  finite  number  of  blades  is  considered. 


The  effect  of  a  finite  number  of  blades  is  to  give  a  further  loss  in  wake  energy  due  to  the 
individual  vortex  sheets  from  each  blade.  Goldstein  worked  out  the  exact  effect  for 
optimized  rotors.  Prandd’s  approximate  correction  factor,  [16],  [17],  [18],  agrees  very  well 
with  Goldstein  for  moderate  climb  rates.  Because  of  this  tip  loss,  for  a  given  thmst,  there  is 
more  induced  flow  than  predicted  by  momentum  theory.  Using  these  principles  an 
approximation  to  the  theoretical  figure  of  merit  for  an  actuator  disk  or  a  lifting  rotor  for  a 
finite  number  of  blades  can  be  obtained  using  Prandd  formuladon.  The  Prandd  k  factor  is 
applied  as  follows. 

dL  =  (inrdr)  p(y  +  v){2v)k 

where 


;  ^  -1 
k  =  — cos 

n 


exp 


IX 


(4.12) 


From  the  above  expression  it  is  seen  that  near  the  root  (r  —>-0)  the  correction  factor  is  close 
to  the  unity,  and  it  decreases  as  the  radial  position  moves  along  the  blade  towards  the  tip. 
This  formula  will  be  used  in  the  results  to  follow. 
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Prandtl’s  result  can  be  used  to  obtain  an  approximate  expression  for  the  figure  of  merit  for 
the  case  of  an  actuator  disk  (lift  perpendicular  to  the  rotor),  in  axial  flow  and  for  the  case  of 
finite  number  of  blades  as  follows. 

do 

Later  in  this  thesis,  we  will  utilize  Eq.  (4.13)  and  compare  with  other  case. 

4.1.3  Betz  Distribution 


When  the  lift  vector  is  tilted  perpendicular  to  the  vortex  sheets,  the  ideal  power  is  no 
longer  attainable.  Thus,  uniform  flow  is  no  longer  the  optimum  condition.  Betz  [18] 
determined  that  the  minimum  power  is  obtained  when  the  induced  flow  at  the  individual 
blades  is  such  that  the  vortex  sheet  remains  along  a  helical  path.  Thus,  the  optimum  inflow 
distribution  is  proportional  to  coscp.  Figure  3.5.  For  an  infinite  number  of  blades,  it  follows 
that  the  pressure  field  must  follow  this  same  shape.  Thus,  let  the  optimum  pressure  at  the 
rotor  be: 


AK  KV  r 

AP  = cos  ©  = - , 

2  2  7771^ 


(4.14) 


where  A  is  a  Lagrange  multiplier.  Then  the  induced  velocity,  2i>  is: 


1  AD  ^ 

IP  = - AP  =  ■ 


21T  4 


(4.15) 


The  thmst  coefficient  is  : 


C !  =  2j^  AP  cos  (prdr  =  KV £  coxxprdr  =  KV - 


-dr 


(4.16) 
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Division  of  both  sides  by  C^,  introduction  of  the  normalized  values  A  =  - 


- ,  and  with  j  =  rlX  we  obtain  an  expression  that  can  be  integrated  to  obtain  a 


close  form  expression  for  the  figure  of  merit  of  the  Betz  distribution. 


AK  2  J 


Jo  1+  / 


(4.17) 


Integration  on  Equation  (4.17)  gives  the  value  of  the  normalized  Lagrange  multiplier  to  be: 


A  =  ± 


^  l-l"ln(  1  +  A 


(4.18) 


To  obtain  the  induced  power  coefficient,  we  must  consider  the  power. 


Cp  =  2^  iSPwrdr  - 


AV  j-i  r 
4  Jor"+;L" 


(4.19) 


Then,  the  normalized  power  coefficient  is: 

-  C  .  . 

Cp^  = - - — jjj,  which  provides  the  final  expression  for  the  normalized  induced  power. 

2  (L;.  /  2^ 


r  J  J 

I 


(4.20) 


—  — 2  ^ 

—  —  —  —  77  77  1 

For  axial  flow,  V  =  X  =ri  +v  =  — h, - h  1  =  . 

2  V  4 


The  result  of  this  closed-form  solution  yields  the  expression  for  the  ideal  figure  of  merit  for 
a  lifting  rotor  with  an  infinite  number  of  blades.  Rearranging  the  previous  equations,  one 


obtains: 
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(4.21) 


F.AV=l-2"lnfl  +  i 

V  ^  y 

It  is  interesting  that  Betz  never  provided  a  solution  for  the  figure  of  merit  of  his 
distribution,  but  provided  all  the  necessary  theoretical  background  to  obtain  this  closed-form 
expression  for  the  efficiency. 

One  can  combine  the  Betz  distribution  with  the  Prandtl  tip-loss  correction  to  obtain  an 
approximate  figure  of  merit  for  tilted  lift  and  tip  losses. 

P -M k  cos"  g)rclr  (4.22) 

This  expression  must  be  determined  numerically. 

4.1.4  Goldstein’s  Distribution 

Prandd  gives  an  approximate  solution  for  the  effect  of  blade  number  on  the  Betz 
optimum  distribution.  The  accuracy  of  this  approximation  is  better  as  the  number  of  blades 
increases  and  as  the  ratio  of  the  tip  speed  to  the  advancing  velocity  becomes  larger  (inverse 
of  the  advance  ratio).  Goldstein  solved  for  the  exact  optimum  circulation  distribution  using 
Betz’s  velocity  distribution  for  an  optimally  efficient  propeller.  He  was  able  to  obtain  the 
exact  solution,  and  that  in  fact  it  was  very  close  to  the  approximate  solution  provided  by 
Prandtl. 
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Figure  4.2  shows  the  results  of  the  two  methods  for  two  and  four-bladed  propellers.  It 
gives  the  circulation  distributions  for  different  velocity-normalized  radial  positions  obtained 
by  Prandtl  and  Goldstein.  Prandtl’s  solution  is  the  dotted  line  and  Goldstein’s  is  the  solid 
curve.  The  circulation  is  higher  for  a  larger  number  of  blades. 


Figure  4.2:  Prantl  and  Goldstein  circulation  distributions  for  a  2  and  a  4-bladed  propeller. 

Tip-speed  ratio  of  5. 


4.2  Finite-State  Methods:  Axial  Flow 


Results  obtained  using  finite-state  methods  cannot  be  shown  without  the  use  of 
previously  introduced  concepts  and  techniques  applied  to  the  case  of  a  Ughtly-loaded  rotor 
that  is  the  subject  of  this  research.  First  of  all,  the  general  optimization  must  be  applied  to 


obtain  the  pressure  coefficients  that  provide  an  optimum  efficient  propeller.  Some  general 
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concepts  have  been  introduced,  but  especially  important  is  to  find  an  expression  for  the 
figure  of  merit  (or  efficiency)  that  we  can  apply  to  all  of  the  cases  in  axial  flow.  This  general 
expression  will  change  according  to  the  assumptions  for  each  of  the  four  cases  studied  in 
axial  flow  (actuator  disk  and  tilted  lift  with  a  finite  or  an  infinite  number  of  blades). 

The  results  will  demonstrate  agreement  with  classical  theories  so  the  finite-state 
methodology  can  be  verified.  For  this  reason,  each  of  the  four  cases  discussed  in  the 
previous  sections  will  be  studied  independendy  and  compared.  The  finite-state  model  will 
be  compared  to:  1)  optimum  pressure  and  inflow  distributions  given  by  momentum  theory 
(for  an  actuator  disk  with  an  infinite  number  of  blades),  2)  figure  of  merit  and  induced 
power  coefficients  given  by  Betz’s  distribution  (for  tilted  lift  and  an  infinite  number  of 
blades),  3)  figure  of  merit  and  induced  power  coefficient  results  given  by  Prandtl’s 
approximation  (for  an  actuator  disk  with  a  finite  number  of  blades),  and  4)  circulation 
distribution  given  by  the  Prandtl  approximation  to  Goldstein  (for  the  case  of  tilted  lift  with  a 
finite  number  of  blades). 


4.2.1  Optimization 


Given  the  general  quadratic  optimization  on  Section  3.3  we  may  now  apply  this  approach 
to  the  case  for  minimum  induced  power  that  we  are  presently  discussing.  For  an  actuator 
disk,  with  infinite  number  of  blades  that  is  lightly  loaded,  we  will  minimize  Cp  for  a  giver  C.p. 
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c. 


y  2ay;  +  y  y  <r, 

n  n  n  t 


subject  to  a  given  Cj. 


(4.23) 


C,.  =  2 


(4.24) 


Physically,  the  coefficients  are  a  Legendre- function  representation  of  V  cos  (p  which 
is  the  loss  of  thrust  due  to  the  tilt  of  the  lift  vector.  The  proof  that  they  are  in  fact  a  fit  is 
further  discussed  on  Appendix  F.  They  are  defined  as: 


cl  =  Jo 

C*  = — J  J  coscpP^  {v^vdv cos{m\i/)d\i/ 

C"  = — f  f  co?,(pP”  {v^vdvdm{my/)dy/ 

where 


for  a  lighdy-loaded  rotor. 


cos^  = 


r  jj.  sin  y/ 


■sJ(r  +  jUsiny/y+Z^ 


(4.25) 

(4.26) 

(4.27) 


(4.28) 


Then,  the  relationship  between  the  coefficients  and  the  function  coscp  (cosine  of  the  tilt 
angle)  becomes: 

vcosy)  =  '^CIp°{v)+'^P”{v)[c^'  cos{my/)  +  C”"  ?.m{my/)]  (4.29) 

n  m^n 
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Notice  that,  for  axial  flow,  the  advance  ratio  ([x)  is  zero.  For  an  ideal  actuator  disk,  the 
non-dimensional  cMmb  rate  (k)  is  arbitrarily  small,  whereas  for  tilted  lift  it  will  have  a  finite 
value. 


From  Fle’s  inflow  Equations  (3.9)  and  (3.10)  for  an  infinite  number  of  blades,  this  is  a 
steady  system.  Furthermore,  all  the  coefficients  associated  with  the  sine  component  are 
zero,  that  is,  fi”  and  r"  are  zero.  Rearranging  the  equation  with  this  in  mind,  it  can  be 
solved  in  matrix  form  for  the  induced  velocity  coefficients,  as  expressed  by: 


F[L-F{a:)=i{r:'} 


2a 


cc. 


1 

2K 


{<) 


Cp=— (y)' 


j  mr 
'^nj 


(4.30) 

(4.31) 

(4.32) 

(4.33) 


Let  the  matrix  L  with  the  m  —  Q  row  partition  multiplied  by  two  will  be  called  L  .  Then, 
the  minimum  induced  power  problem  can  be  formulated  as  the  optimization  of  a  functional 
J,  shown  below,  subject  to  the  constraint  expressed  by  equation  (4.24).  The  optimization  is 
carried  out  in  to  yield  the  optimum  value  for  the  pressure  coefficients. 


/  =  -  — 
2  K 


jjn 


(4.34) 
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subject  to  the  condition 


1^,  I  =  Q 


This  leads  to 


d 


1  1 

2  K 


(r*V  iC”] 


=  0 


2  " 


I-;  !*;)-a{*;)  {c; 


=  0 


1 

K 


j  nn 
'^nj 


1 

+  - 
2 


{C)={C}A 


V  )  optimal 


j  rm 
^nj 


1 

+  - 
2 


J  rm 
^nj 


{c:)a=k[l,J'{c:)a 


where  A  is  chosen  to  match  the  desired  C,  in  Equation  (4.35). 


(4.35) 

(4.36) 

(4.37) 

(4.38) 

(4.39) 


Equation  (4.39)  is  the  solution  for  this  optimization  problem,  which  will  yield  the 
minimum  induced  power  for  a  lighdy  loaded  actuator  disk  with  infinite  number  of  blades. 

In  these  equations  A  is  the  Lagrange  multiplier,  that  is  chosen  to  give:  r®  =  . 


4.2.2  Figure  of  Merit 


The  general  solution  for  the  pressure  coefficients  given  in  Section  4.2.1  can  be  applied  to 
different  cases.  It  is  the  purpose  of  this  section  to  show  results  for  axial  flow,  but  these 
coefficients  can  be  also  used  to  obtain  pressure,  circulation  and  inflow  velocity  for  a  variety 
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of  flows,  including  edgewise  flow  (j  —  90°).  For  axial  flow,  j  —  0°,  the  elements  in  L  are 
zero  except  when  r  —  m. 


For  an  infinite  number  of  blades,  the  equations  of  motion  reduce  to 


(4.40) 


The  pressure  and  induced  velocity  coefficients  can  then  be  obtained  and  introduced  into  the 
equation  for  the  power  and  thrust  coefficients. 


{<}  =  2v[LJ'{a:}  (4.41) 

=  (“■«) 


The  power  coefficient  becomes 


\2a^' 

.4  .  irl 

}  « 

II 

1  fl 

f 

/ 

\ _ 

/ - 

(4.43) 


For  this  optimization,  we  find  the  power  subject  to  a  thrust  constrained,  expressed  in 
Equation  (4.44). 

C,=2{c}'{c)  (4.44) 


Introducing  the  induced  velocity  coefficients  in  Equation  (4.45)  into  the  equation  for  the 
power  (Equation  (4.43))  provides  the  power  to  minimize. 

1 


I  2a 
I  a 


0 


2K 


E 


{<} 


(4.45) 
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Cp— {C}'[C]{C}=— {cr[£,]!' 


(4.46) 


For  minimum  power,  without  the  tilt  on  the  Uft,  the  optimum  pressure  and  influence 


coefficients  are  =  l/Vs  and  r”  =  . 


The  functional  for  the  optimization  becomes 


(4.47) 


Performing  the  optimization  (dj  —  0)  the  distribution  that  provides  optimum  power  is 


{c}=c[l,J‘{c:}. 


(4.48) 


For  the  case 


ofno  tilt,  r”  =(73/2) Q  =K  L  '(i/V3)a. 


The  first  value  of  the 


influence  coefficient  matrix  obtained  using  the  closed-form  expressions  in  Appendix  D  is 


Ls  =3/2  and  Ls  =3/4. 


The  power  coefficient  and  the  Lagrange  multiplier  are 


(4.49) 


CjjlV 


(4.50) 
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Introducing  Eqs.  (4.49)  and  (4.50)  into  the  equation  for  the  power  coefficient  yields  the 
optimum  power  coefficient  in  Equation  (4.51). 


L 


L 


L 


\  «  )  OTX 


2K 


v2Ky 


{c)'TrT'{c) 


r  r  ^  7' 

-  —  n 

I-E  ^  1 

1 - 

I _ 

^  (4.51) 


Multiplying,  the  final  form  of  the  optimum  power  coefficient  is 

C,," 


8V 


{c} 


(4.52) 


Applying  the  same  concepts  to  the  case  of  axial  flow  with  no  tilt  on  the  lift,  where 
|C*}  =  |l/  Vs  0  0  ...|  the  minimum  power  coefficient  is  obtained. 


c  ' 

'^7 

'  1 

L 

-1  1 

8V 

_V3 

s  1 

”  V3_ 

'  c  ^ 

_  I 


2K 


(4.53) 


Introducing  these  power  coefficients  into  the  definition  of  figure  of  merit  demonstrates  the 
steps  taken  to  come  to  the  final  expression  used  throughout  this  dissertation. 


PM  =  =  4{c; }'"  [ls\  {c  } 

E  p 


(4.54) 


4.2.3  Infinite  Number  of  Blades 


For  the  case  of  an  actuator  disk.  Equation  (4.28)  reduces  to  the  unity  (X  —  0)  and  the 
coefficients  in  equations  (4.25)  through  (4.27)  reduce  to: 
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dv  =  ^ 

s 

and  all  the  others  become: 

C=o 

Therefore,  the  vector  is  the  vector  {l/Vs  0  0  ...  with  as  many  elements  as  the 

number  of  terms  that  correspond  to  the  harmonics  studied  in  the  problem 

Momentum  theory  predicts  that  the  minimum  induced  power  for  an  actuator  disk  in  axial 
flow  will  be  achieved  by  constant  pressure  and  constant  inflow  distributions.  Results  using 
the  finite-state  method  show  agreement  with  these  predictions.  Figures  4.3  and  4.4  show  the 
constant  profile  for  the  pressure  and  the  inflow  respectively  with  a  reduced  amount  of  terms 
in  the  Fourier  series.  According  to  blade-element  theory,  the  lift  distribution  corresponding 
to  this  pressure  and  velocity  should  be  linear.  Finite-state  methods  agree  with  the  predicted 
results,  as  it  is  shown  in  Figure  4.5. 

Figures  4.6  through  4.17  show  the  pressure  and  induced  velocity  distributions  for  all 
azimuthal  locations  on  the  rotor  for  increasing  values  of  cUmb  and  a  thrust  coefficient  of  1.0. 
The  pressure  and  induced  velocity  for  axial  flow  and  the  case  for  no  tilt  (u=0)  are  constant 
for  all  angles  around  the  rotor,  as  predicted  from  classical  theories.  For  this  case  only  (of  no 
climb)  the  rotor  is  in  hover,  therefore  the  assumption  that  it  is  lightly-loaded  no  longer 
holds.  Both  pressure  and  velocity  are  directly  related,  therefore  it  is  predicted  that  they 
behave  similarly  for  increasing  values  of  the  tilt.  It  is  seen  that  an  adverse  inflow  region 
appears  at  the  center  of  the  rotor  and  that  this  region  augments  with  increasing  tilt.  As  the 


c‘=\'Sv- 
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pressure  and  inflow  drop  the  flgure  of  merit  decreases  (i.e.,  the  efficiency).  This  decreased 
efficiency  will  be  studied  more  in  detail  when  discussing  the  forward  flight  results. 

One  way  to  make  sure  that  the  theory  in  fact  provides  results  that  were  expected,  apart 
from  the  fact  that  for  the  simplest  case  both  pressure  and  inflow  are  constant,  is  to  see  how 
the  curvature  of  the  distributions  at  the  center  of  the  rotor  behaves.  As  the  depressed  region 
at  the  root  increases,  previous  theories  predict  that  the  pressure  and  inflow  should  increase, 

close  to  the  root,  as  juj ^Jl  +  jU^  .  Plots  for  any  individual  a2imuthal  location  and  a  given  tilt 
show  that  these  results  follow  the  predicted  shape  (See  Refs.  [19]  and  [10]). 

Of  importance  is  to  notice,  in  all  of  these  figures,  that  the  air  is  coming  into  the  rotor  at 
an  angle.  The  rotor  advancing  side  is  where  rcostf)  =  -1.  Therefore  the  aft,  where  ([j  =  0 
corresponds  with  rcos^  =  1.  There  is  a  symmetry  for  axial  flow,  but  this  symmetry  will  no 
longer  hold  for  skewed  flow  and  therefore  the  direction  to  which  the  rotor  is  advancing  will 
become  important. 
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Figure  4.3:  Pressure  profile  that  provides  minimum  induced  power  for  an  actuator  disk  in  axial  flow 

with  an  infinite  number  of  blades. 
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Figure  4.4:  Velocity  profile  that  provides  minimum  induced  power  for  an  actuator  disk  in  axial  flow 

with  an  infinite  number  of  blades. 


Figure  4.5:  Lift  distribution  for  optimum  flow  for  an  actuator  disk  with  an  infinite  number  of  blades. 
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Pressure  distribution  for  u  =  0 


r*cos(i;;) 

Figure  4.6:  Pressure  Distribution, )(  =  0,  u  =  0.0,  FM=1.0. 


Pressure  distribution  for  x=0.  u  =  0.2 


r*cos(v);) 


Figure  4.7:  Pressure  Distribution,  ^  =  0,  u  =  0.2,  FM=0.87. 
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Pressure  distribution  for  x=0.  u  =  0.4 


r*cos(nr) 

Figure  4.8:  Pressure  Distribution,  ^  =  0,  u  =  0.4,  FM=0.69. 


Pressure  distribution  for  x=0,  u  =  0.6 


r^cooft);) 

Figure  4.9:  Pressure  Distribution,  x  =  0,  u  =  0.6,  FM=0.52. 
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Pressure  distribution  for  x=0,  u  =  0.8 


r*cos(v) 

Figure  4.10:  Pressure  Distribution,  x  =  0,  u  =  0.8,  FM=0.40. 


Pressure  distribution  for  x=0,  u  =  1.0 


r^cDsfi};) 

Figure  4.11:  Pressure  Distribution,  ^  =  0,  u  =  1.0,  FM=0.31. 
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V*lnduced  velocity  distribution  for  x=0,  u  =  0 


r*cos(Hr) 


Figure  4.12:  Induced  Velocity  Distribution,  x  =  0,  u  =  0.0,  FM=1.0. 


Induced  velocity  distribution  for  x=0.  u  =  0-2 
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Figure  4.13:  Induced  Velocity  Distribution,  x  =  0,  u  =  0.2,  FM=0.81. 
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Induced  velocity  distribution  for  x=0.  u  =  0  4 


r^cosCv);) 


Figure  4.14:  Induced  Velocity  Distribution,  y  =  0,  u  =  0.4,  FM=0.69. 


Induced  velocity  distribution  for  x=0,  u  =  0.6 


r*cos(n;) 


Figure  4.15:  Induced  Velocity  Distribution,  y  =  0,  u  =  0.6,  FM=0.52. 
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Induced  velocity  distribution  for  x=CI.  u  =  0-8 


r’cosfv) 

Figure  4.16:  Induced  Velocity  Distribution,  x  =  0,  u  =  0.8,  FM=0.40. 


Induced  velocity  distribution  for  x=0.  u=10 


r*cos(i|;) 


Figure  4.17:  Induced  Velocity  Distribution,  x  =  0,  u  =  1.0,  FM=0.31. 
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Since  the  above  shows  that  a  formal  optimization  with  the  dynamic  wake  model  gives  the 
Glauert  result  for  minimum  power,  it  seems  that  it  would  be  useful  to  consider  some  closed- 


form  results  under  the  Glauert  hypothesis.  From  momentum  theory  one  can  show  that 


Cj  =  lirj  -1-  v)v 

(4.55) 

Cp  =2{'q  +  vfv  =  {'q+v)C-p 

(4.56) 

Because  we  optimize  for  a  given  C.J.,  it  is  very  convenient  to  normalize  all  velocities  on 
induced  flow  in  hover.  Thus,  rj  =7]  !  ,  V  =  v !  ^JCj.  /  2  ,  A  =  rj  +v  .  It  follows  that 

the  proper  normalization  of  induced  power  is 


^  _^/2G 

^  D 


C  ' 

'^7' 


3/2 


(4.57) 


The  result  is  a  normalized  set  of  inflow  equations.  The  thrust  equation  becomes: 

\  =  {ff  +  v)v  (4.58) 

which  can  be  solved  for  normalized  or  flow  due  to  a  normalized  cMmb  rate.  That  value  can 
then  be  used  to  determine  the  normalized  induced  power  for  an  ideal  actuator-disk  rotor. 


Cp  =Cp-ri  = 


f  —2  ^ 

1/2 

2 

l4  y 

(4.59) 


The  ideal  induced  power  ranges  from  a  normalized  value  of  unity  in  hover  =  O)  and  then 
decreases  with  climb  rate  as  1  /  ?/  .  This  is  the  Glauert  result.  We  will  use  this  ideal  value  to 
compare  minimum  power  settings  for  various  rotors.  We  will  define  a  generalized  figure  of 
merit  which  is  the  ideal  power,  equation  (4.33),  divided  by  the  actual  induced  power. 
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The  question  remains  if  the  use  of  finite-state  methods  will  suffice  to  obtain  the  Betz 
distribution  for  a  lifting  rotor  with  an  infinite  number  of  blades.  To  verify  this,  the  figure  of 
merit  is  found  by  finite-state  methods  using  the  formulation  described  in  the  optimization 
section.  These  coefficients  represent  the  general  solution.  To  customize  them  to  the 
present  special  case,  some  modifications  were  performed.  Since  the  present  cases  are  for 


axial  flow,  the  matrix 


—  defined  by  equation  (3.12)  —  simplifies  to  a  diagonal  in  terms 


of  r”  .  It  should  be  noted  that  these  matrices  are  identical  to  in  Ref.  4.  The  closed 

form  expression  for  the  matrix  can  be  seen  in  Appendix  D.  For  an  infinite  number  of 
blades  in  axial  flow,  C”  =  0  except  when  =  0  so  that  only  enters  the  optimization. 


comes  from  the  following  integral  over  the  wake  skew  angle. 


C 


0 

n 


(k) 


where  cos^=  ,  and  rdr=-vdv.  With  these  influence  coefficients  the  thrust  and 

power  coefficients  become: 

Cr=2  Y.  C.t":  (4.60) 

C,=riC,+  X  (4.61) 


For  this  optimization,  again  the  power  coefficient  is  minimized  for  constant  thrust,  and  the 
functional  /becomes: 

i  =  ''C,  +  I  Z  C.t:  (4.62) 
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where  A  is  the  Lagrange  multiplier.  Performing  the  optimization  (8/  —  0),  the  optimal 
pressure  coefficients  for  this  particular  case  are: 


{""U„=[AT'L.}AF  (4.63) 

Introducing  the  above  changes  to  the  general  optimization  formulation,  the  figure  of  merit 
using  finite-state  methods  is: 

PM., =  2{C,}'  [4]^'  {C.}  (4.64) 


Figure  4.18  shows  the  comparison  for  the  figure  of  merit  from  the  finite-state  method  as 
compared  to  the  Betz  formula.  It  is  seen  that  the  finite-state  method  agrees  satisfactorily 
with  Betz  result.  The  difference  between  them  can  be  reduced  by  addition  of  more  terms  to 

I^Ll”  ^  and  .  However,  the  present  approximation,  which  uses  twenty  terms  is  thought 

to  be  close  enough  so  that  the  dynamic  inflow  model  is  verified. 

What  is  most  important  about  Figure  4.18  is  the  large  drop  in  figure  of  merit  with  climb  rate, 
even  for  an  optimized  rotor.  The  drop  is  due  purely  to  the  effects  of  tilted  lift  and  swirl 
velocity.  It  may  well  be  that  the  deficiency  in  rotor  efficiency  in  forward  flight  is  due  to  a 
similar  phenomenon. 

Figure  4.19  gives  figure  of  merit  as  a  function  of  cMmb  rate  r\  rather  than  X.  Since  X  is 
determined  by  the  total  flow  through  the  rotor,  the  Figure  of  Merit  thus  becomes  a 
function  both  of  climb  rate  r\  and  thrust  coefficient  C^..  For  larger  cUmb  rate,  the  effect  of 
C.J.  is  diminished. 
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Figure  4.20  presents  the  induced  power  coefficient  as  a  function  of  climb  rate  for  a  range 


of  C.J.  values 
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Figure  4.20:  Induced  power  coefficient  (helicopter  convention)  versus  climb  rate  (infinite  number  of 

blades) 
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An  important  characteristic  seen  in  Figure  4.20  is  a  “bucket”  in  the  induced  power  for 
each  thmst  coefficient  at  a  given  climb  rate.  This  is  due  to  the  fact  that  ideal  power 
decreases  with  r\  whereas  the  figure  of  merit  also  decreases  with  r].  Thus,  there  is  an 
optimum  cHmb  rate.  The  lowest  curve  is  for  C.j.  =  0  and  is  equal  to  This  ideal 

minimum  power  monotonically  decreases  with  rj ,  so  the  “bucket”  is  not  present,  and  the 
induced  power  coefficient  does  not  increase  for  high  cMmb  rates. 


4.2.4  Finite  Number  of  Blades 


Because  the  Prandd  correction  as  applied  to  the  Betz  distribution  agrees  so  closely  with 
Goldstein,  that  it  makes  sense  to  do  some  calculations  with  the  Prandd  factor  to  determine 
the  magnitude  of  the  effect  of  number  of  blades  on  figure  of  merit.  Thus,  equadon  (4.12) 
can  be  used  for  the  figure  of  merit  computadons. 

Figure  4.21  shows  the  effect  of  dp  loss  (as  determined  from  Prandd’s  /^-factor)  on  figure 
of  merit.  The  top  curve,  for  the  case  cp  =  0,  is  the  effect  for  an  actuator  disk  with  a  finite 
number  of  blades  for  a  four-bladed,  Ughdy  loaded  rotor.  The  middle  two,  coincident  curves 
are  the  figure  of  merit  for  a  rotor  with  tilted  lift  but  infinite  number  of  blades.  The  lowest 
curve  is  for  dlted  lift  and  finite  number  of  blades  (i.e.,  the  Goldstein  solution).  One  can  see 
that  blade  number  is  also  an  important  factor  in  the  loss  of  ideal  induced  power. 
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of  Merit 


Figure  of  Merit  Comparison  for  any 


Figure  4.21:  Effect  of  tip  loss  on  figure  of  merit,  lightly-loaded  rotor. 
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We  now  wish  to  see  if  the  finite-state  methodology  can  give  the  correct  optimum 
distribution  and  figure  of  merit  as  Goldstein  (i.e.,  as  the  Prandd-corrected  Betz).  Makinen 
showed  that  the  inflow  model  can  indeed  match  Goldstein  provided  that  a  correction  is 
applied  for  the  swirl  kinetic  energy.  Thus,  the  added  energy  is  added  to  the  mass  matrix;  and 
the  resultant  induced  flow  is  assumed  parallel  to  the  tilted  lift  vectors. 

To  be  precise,  the  apparent  mass  matrix  (diagonal),  must  be  replaced  to  include  the 

effect  of  the  wake  swirl.  There  are  different  swirl  corrections  that  can  be  applied,  but  from 
Ref.  10  the  following  correction  gives  the  best  results. 

r  /  -  1 

[k]^\4k]  f]+»  ^  [w-G-dT  \'l^]  <'*■“) 

L  J 

where  k  =2.2,  ^  is  the  number  of  blades,  and  X  is  the  total  inflow.  It  should  be  noted  that 
for  an  actuator  disk  (no  lift  tilt),  k  is  set  to  zero. 

Performing  the  optimization  for  this  case,  and  using  finite-state  methods,  the  figure  of 
merit  is  obtained.  The  vector  |G*|  remains  the  same  as  for  an  infinite  number  of  blades. 

There  is  no  change  in  the  values  because  physically  these  coefficients  are  a  fit  of  the  function 
V  cos  (p ,  and  for  an  actuator  disk  cos  (p  =  \. 

F-M.f,._,,.=2{c:Y[Ly{c:]  (4.66) 

Once  the  theory  has  been  verified,  some  useful  plots  of  induced  power  for  different 
numbers  of  blades  at  various  climb  rates  can  be  obtained,  as  it  is  shown  in  Figure  4.22.  The 
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importance  of  this  graph  is  that  the  effect  of  finite  number  of  blades  on  the  induced  power 
can  be  noticed.  It  is  seen  that  induced  power  increases  for  a  decreasing  number  of  blades. 
It  is  an  expected  result,  as  the  ideal  induced  power  exists  for  an  infinite  number  of  blades 
(for  Prandd  is  /fe  =  1).  The  profile  of  the  curves  is  similar  to  the  one  observed  for  infinite 
numbers  of  blades  at  different  thrust  coefficients.  The  “bucket”  effect  is  present  here  also, 
and  the  general  profile  is  maintained.  Thus,  the  effect  of  these  differences  for  finite  number 
of  blades  affecting  the  induced  power  is  as  less  critical  as  the  increase  in  induced  power  due 
to  lift  hit. 
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Finite-State  methods  should  agree  with  the  theory  developed  by  Goldstein  [1 9]  for  every 
flight  condition  in  axial  flow.  There  is  no  closed-form  solution  or  expression  that  Goldstein 
developed  for  the  theoretical  figure  of  merit  for  a  lifting  rotor  with  a  finite  number  of  blades. 
However,  Makinen  [10], [11]  was  successful  in  the  further  development  and  application  of 
the  finite-state  method  to  obtain  circulation  for  a  given  induced  velocity.  These  circulation 
results  are  in  agreement  with  Goldstein’s  circulation  for  an  optimal  propeller,  as  it  is  shown 
in  Figures  4.23  and  4.24.  The  fact  that  the  application  of  finite-state  methods  provides  an 
accurate  optimal  circulation  results  in  the  confidence  that  the  calculations  of  figure  of  merit 
for  this  special  case  will  also  be  accurate. 

Figures  4.23  and  4.24  show  the  circulation  at  any  radial  position  of  the  blade  using 
Prandtl’s  approximation,  Goldstein’s  optimal  circulation,  and  Makinen’s  results  with  the 
swirl  velocity  corrections  made  to  the  apparent  mass  matrix  in  equation  (4.48).  Figure  4.23  is 
for  a  //q  =  5  (A  =  0.20)  and  Figure  4.24  is  for  =  20  (A  =  0.05).  It  is  noticed  that  Prandt  and 
Goldstein’s  circulations  give  results  that  are  very  close  to  each  other.  Since  there  is  such 
close  agreement  in  both  approaches,  and  there  is  a  figure  of  merit  expression  for  Prandt’s 
approximation,  the  finite-state  approach  could  be  comparable  to  Prandd’s  approximation.  It 
is  not  surprising  that  the  quadratic  optimization  with  the  dynamic  wake  model  gives  the 
correct  figure  of  merit  due  to  both  lift  tilt  and  finite  number  of  blades.  Figures  4.23  and  4.24 
(from  Ref.  8)  show  that  the  dynamic  wake  model  (with  swirl  correction)  gives  the  correct 
inboard  (swirl)  and  outboard  (tip  loss)  velocities. 
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Circulation 


Figure  4.23:  Circulation  at  any  blade  radial  location  for  Prandtl,  Goldstein,  and  using  Finite-State 

methods.  Plot  obtained  from  [10], [11]. 
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Circulation 


Figure  4.24:  Circulation  at  any  blade  radial  location  for  Prandtl,  Goldstein,  and  using  Finite-State 

methods.  Plot  obtained  from  [10], [11]. 
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4.3  Finite-State  Methods:  Skewed  Flow 


Since  the  method  has  been  validated,  the  theory  can  be  applied  to  the  same  cases  for 
skewed  flow.  Most  of  the  approach  for  skewed  flow  is  similar  to  that  for  axial  flow. 


4.3.1.  Infinite  Number  of  Blades 


The  case  for  an  actuator  disk  with  an  infinite  number  of  blades  will  be  revisited. 
However  some  changes  will  be  done  for  forward  flight.  For  this  case,  momentum  theory 
predicts  that  the  optimum  induced  power  is  obtained  for  a  constant  pressure  distribution 
(similarly  as  to  what  happened  for  axial  flow)  but  the  induced  velocity  profile  will  no  longer 
be  constant.  For  an  actuator  disk  with  an  infinite  number  of  blades,  we  have  already  applied 
the  finite-state  model  and  verified  that  it  gives  the  Glauert  solution  of  uniform  pressure. 
However,  to  go  on  to  the  other  cases,  all  harmonics  (and  their  periodic  coupling)  will  need 
to  be  included. 

The  rest  of  the  special  cases,  for  an  actuator  disk  with  a  finite  number  of  blades  and  for 
the  two  cases  for  a  lifting  rotor,  will  provide  results  never  obtained  before.  The  results  will 
hopefully  provide  the  conclusion  as  to  why  the  experimental  minimum  induced  power  for  a 
helicopter  is  orders  of  magnitude  greater  to  what  theory  predicts  should  be.  These  results 
could  allow  determining  what  changes,  if  any,  should  be  introduced  in  the  rotor  to  reduce 
the  minimum  induced  power. 
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The  formulation  for  the  figure  of  merit  in  forward  flight  will  remain  similar  to  the  general 
figure  of  merit  shown  by  equation  (4.54): 

However,  the  coefficients  and  the  L-matrix  will  be  different  than  the  ones  obtained 
previously,  and  also  different  for  each  of  the  four  cases. 

The  main  difference  for  skewed  flow  is  that  when  calculating  the  cosine  of  the  inflow  angle 

the  advance  ratio,  g,  must  be  considered.  Equation  (4.28)  again  is: 

r  +  ju  sin  ur 
cos(p=  .  ^ 

+  jUsini//^  + 

where  r  is  the  radial  position  along  the  blade  and  y/  is  the  angle  at  which  the  rotating  blade 
is  with  respect  to  the  aft  position  of  the  rotor. 

For  skewed  flow,  the  total  inflow  also  changes.  The  total  inflow  for  axial  flow  was  defined 
before  as: 

A  =ri  +v  =  —  + 

2 

and  it  was  derived  from  momentum  theory  for  a  uniform  induced  flow  distribution.  In 
forward  flight,  the  total  inflow  becomes: 

A  =fj  +  v  (4-67) 

where  the  normalHed  inflow  is  the  solution  of  equation  (4.53)  for  given  normalized  climb 
rate  and  advance  ratio. 
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These  changes  will  affect  the  optimum  coefficients,  but  the  L-matrix  will  also  be  altered 
because  the  skew  angle  is  no  longer  zero,  and  so  there  are  more  harmonics  than  the  m  —  ^ 
for  axial  flow.  The  expression  for  this  matrix  will  be  obtained  using  He’s  formulation  (Eqs. 
(3.9)  and  (3.10)). 

Different  cases  have  been  considered  for  a  rotor  with  low  lift  and  with  the  tilt  angle 
defined  above.  Table  4.1  shows  all  of  the  cases,  including  the  axial  flow  cases  shown  in 
Section  4.2.3  (Figures  4.6  through  4.17).  The  Table  provides,  for  each  case  the  value  of  the 
skew  angle  used,  the  advance  ratio,  the  total  inflow  and,  based  on  these  parameters,  the  type 
of  case  that  the  results  show.  None  of  the  cases  considers  the  swirl  correction,  because  it  is 
small  in  forward  flight.  Axial  flow  solutions  (Section  4.2.3)  are  cases  1  through  7.  Case 
number  1  is  for  no  tilt  and  cases  2  through  7  are  for  a  climbing  actuator  disk  but  a  disk  that 
is  not  traveling  forward  (i.e.,  no  advance  ratio).  Cases  7  through  28  are  for  various 
conditions  in  forward  flight.  Of  special  attention  are  cases  8,  15  and  22  that  show  pure 
skewed  flow,  with  no  hit  or  advance.  Cases  23  through  28  show  skewed  flow  advancing 
without  climb. 

For  cases  with  low  lift,  the  mass  flow  parameter,  V  becomes  u.  The  theory  developed 
enables  anyone  to  obtain  results  for  a  case  with  lift  just  by  introducing  the  modified  values 
that  affect  the  mass  flow  parameter  and  the  total  inflow,  X.  The  equations  used  to  obtain  the 
values  of  advance  ratio  and  total  inflow  for  all  of  the  cases  in  Table  4.1  are  shown  in 
Equations  (4.69)  through  (4.72). 
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+{A  +  v)A. 

(4.69) 

«  =  90°  -  j 

(4.70) 

ju  =  t>cos(«) 

(4.71) 

Z  =  usin[a)  +  v 

(4.72) 

where  v  —  0  for  low-lift  cases,  x  is  the  skew  angle  and  u  is  the  total  inflow. 

The  cases  on  Table  4.1  are  applied  to  figure  of  merit  calculations  and  to  obtain  pressure 
and  induced  velocity  distributions  for  all  of  the  points  around  the  rotor. 

Figure  4.25  shows  the  figure  of  merit  for  all  skew  angles.  Figure  4.26  shows  the  figure  of 
merit  change  with  increasing  cHmb  and  advance  ratio. 
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Table  4.1:  Cases  for  FM,  Pressure  and  Inflow  Figures. 


CASE: 

o 

;r(DE 

G) 

1 

TYPE: 

1 

0.0 

0 

0.0000 

0.0000 

Actuator  disk,  axial  flow. 

No  flit 

2 

0.1 

0 

0.0000 

0.1000 

3 

0.2 

0 

0.0000 

0.2000 

4 

0.4 

0 

0.0000 

0.4000 

5 

0.6 

0 

0.0000 

0.6000 

6 

0.8 

0 

0.0000 

0.8000 

7 

1.0 

0 

0.0000 

1.0000 

8 

0.0 

30 

0.0000 

0.0000 

Actuator  disk,  skewed  flow 

No  tilt,  no  advance 

9 

0.1 

30 

0.0500 

0.0866 

Actuator  disk,  skewed  flow 

10 

0.2 

30 

0.1000 

0.1732 

11 

0.4 

30 

0.2000 

0.3464 

12 

0.6 

30 

0.3000 

0.5196 

13 

0.8 

30 

0.4000 

0.6928 

14 

1.0 

30 

0.5000 

0.8660 

15 

0.0 

60 

0.0000 

0.0000 

Actuator  disk,  skewed  flow 

No  tilt,  no  advance 

16 

0.1 

60 

0.0866 

0.0500 

Actuator  disk,  skewed  flow 

17 

0.2 

60 

0.1732 

0.1000 

18 

0.4 

60 

0.3464 

0.2000 

19 

0.6 

60 

0.5196 

0.3000 

20 

0.8 

60 

0.6928 

0.4000 

21 

1.0 

60 

0.8660 

0.5000 

22 

0.0 

90 

0.0000 

0.0000 

Actuator  disk,  skewed  flow.  No  tilt,  no 
advance 

23 

0.1 

90 

0.1000 

0.0000 

Actuator  disk,  skewed  flow 

No  tilt 

24 

0.2 

90 

0.2000 

0.0000 

25 

0.4 

90 

0.4000 

0.0000 

26 

0.6 

90 

0.6000 

0.0000 

27 

0.8 

90 

0.8000 

0.0000 

28 

1.0 

90 

1.0000 

0.0000 
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Figure  4.26:  Figure  of  Merit  change  with  Total  Inflow.  20  harmonics. 


Figure  4.28:  Figure  of  Merit  change  with  Total  Inflow.  1  harmonics. 
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Figure  4.29:  Figure  of  Merit  change  with  Skew  Angle.  20  harmonics,  constrained  inflow  states 


Figure  4.30:  Figure  of  Merit  change  with  Total  Inflow.  20  harmonics,  constrained  inflow  states 
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Figure  4.26  shows  the  figure  of  merit  change  for  different  skew  angles.  For  the  case  of 
edgewise  flow  (^90°),  the  figure  of  merit  does  not  change  with  u  because  there  is  no  tilt,  X. 
To  understand  what  is  happening  with  this  curve,  we  can  see  the  definition  of  the  cosine  of 
tilt  angle,  coscp  in  Equation  (4.28).  For  any  advance  ratio,  if  ^=0,  coscp  equals  1.  Therefore, 
the  influence  coefficient  vectors  for  all  of  these  cases  are  the  same,  which  make  the  figure  of 
merit  equal  for  all  cases.  Physically,  the  fact  that  the  figure  of  merit  does  not  change  for 
edgewise  flow  is  predictable  since  all  of  the  incoming  flow  is  oriented  on  the  same  plane  of 
the  rotor,  so  there  are  no  losses  on  the  perpendicular  direction. 

The  tendency  for  the  figure  of  merit  for  the  rest  of  the  skew  angles  in  Figure  4.26  is  given  by 
a  combination  of  tilt  and  advance.  For  the  case  of  axial  flow  Cl^O”),  the  loss  of  efficiency  is 
the  greatest.  For  axial  flow,  the  advance  ratio  is  zero,  therefore  all  of  the  reduction  in  figure 
of  merit  is  due  to  tilt. 

Several  test  points  are  taken  and  compared  to  previous  results  and  validate  the  method  to 
obtain  the  curves  in  Figure  4.26.  Figure  4.21  is  used  to  check  for  the  case  of  axial  flow 
(lowest-ordinate  curve  in  Figure  4.26).  Figure  4.21  provides  the  figure  of  merit  change  with 
total  inflow  for  different  cases.  The  curves  for  Betz  and  Finite-state  methods  for  infinite 
number  of  blades  are  used  to  check  the  curve  for  2^0°  of  Figure  4.26.  For  \)=1,  ^=1.  For 
example,  the  value  of  figure  of  merit  for  this  point  in  Figure  4.26  is  0.3.  On  the  curve  for 
Betz’s  figure  of  merit  in  Figure  4.21,  for  a  ^=1,  the  value  is  also  0.3.  This  result  proves  that 
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the  curves  for  Figure  4.26  are  accurate  and  depict  the  change  in  figure  of  merit  for  different 
conditions. 

Figures  4.27  and  4.28  show  the  same  as  the  above  but  for  the  case  of  one  harmonic.  The 
number  of  terms  is  drastically  reduced  from  231  to  3  (due  to  the  harmonics  and  the 
corresponding  inflow  states).  A  table  to  determine  the  inflow  states  for  any  umber  of 
harmonics  is  shown  in  Appendix  E.  To  understand  the  effect  of  introducing  a  larger 
number  of  harmonics,  we  take  the  rotor  and  divide  it  radially  and  azimuthaUy  into  sections. 
The  more  harmonics,  the  more  azimuthal  stations  and  the  more  inflow  states,  the  more 
radial  partitions.  Because  of  the  reduced  number  of  partitions,  it  is  not  surprising  that  the 
accuracy  of  the  calculations  is  greatly  reduced.  Flence  the  reduction  of  the  figure  of  merit 
reading  for  all  of  the  cases.  Figures  4.27  and  4.28  demonstrate  the  importance  of  the  use  of 
more  harmonics  and  their  corresponding  inflow  states  in  the  calculations.  A  similar 
conclusion  was  drawn  for  the  case  of  axial  flow  prior  to  making  the  necessary  calculations 
for  forward  flight.  The  determination  of  the  appropriate  number  of  harmonics  is  discussed 
more  in  detail  on  Appendix  E.  If  only  one  harmonic  is  considered,  the  optimization  would 
be  performed  taking  into  account  that  only  one  characteristic  can  be  varied.  For  example,  if 
only  collective  pitch  could  be  controlled,  then  once  the  optimization  is  performed,  the 
efficiency  values  on  Figures  4.27  and  4.28  would  be  the  result  of  this  only  variation.  It  is 
seen  that  the  figure  of  merit  decreases  considerably  if  only  one  characteristic  could  be  varied. 

Finally,  Figures  4.29  and  4.30  are  calculated  to  check  the  impact  on  figure  of  merit  when 
the  inflow  states  are  constrained.  For  these  cases  the  number  of  harmonics  is  chosen  to  give 
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accurate  results,  but  the  number  of  inflow  states  is  constrained.  Physically,  the  pressure 
distribution  (through  the  pressure  coefficients,  r”'s)  is  constrained  to  only  allow  non-zero 

values  for  the  zeroth  and  the  first  harmonic  {m—^  and  m—\  for  sine  and  cosine  only).  The 
impact  on  figure  of  merit  is  not  perceptible,  which  means  that  we  could  disregard  the  rest  of 
the  partitions  as  they  do  not  contribute  on  figure  of  merit.  However,  the  repercussion  on 
pressure  and  induced  flow  would  be  greater,  since  some  of  the  areas  where  the  distributions 
have  abrupt  changes  would  not  be  depicted  as  accurately. 


Figures  4.31  through  4.52  show  the  pressure  and  velocity  distributions  for  all  of  the  cases 
in  Table  4.1  that  correspond  to  skewed  flow.  Note  that  the  pressure  distribution  in  Figure 
4.55  and  the  induced  velocity  distribution  shown  in  Figure  4.56  correspond  to  edgewise  flow 
for  all  values  of  tilt  and  advance.  For  the  induced  velocity  plots,  the  shape  for  all  will  be  the 
same,  but  the  magnitude  will  change  proportional  to  t).  When  the  skew  angle  is  90°,  then 
A=  0.  Going  back  to  the  equation  for  the  cosine  of  the  tilt  angle,  if  the  cUmb  is  zero,  then 
cos  ^3  =  1 ,  independendy  of  the  advance  ratio  /u .  Therefore,  the  pressure  and  velocity 
distributions  do  not  change. 

It  has  already  been  explained  how  the  inflow  and  pressure  distributions  change  for  axial 
flow.  In  general,  for  increasing  tilt  and  advance,  the  distributions  at  the  root  of  the  rotor 
decreases  with  increasing  tilt.  A  similar  pattern  appears  for  all  of  the  cases  in  forward  flight 
but  with  a  main  difference.  In  forward  flight,  as  t)  increases,  and  the  skew  angle  is  increased. 
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the  depression  at  the  root  shifts  to  the  reversed-flow  region  and  becomes  downwash 
(negative). 

Figure  4.27  shows  the  pressure  distribution  around  the  rotor.  As  predicted  by 
momentum  theory,  this  optimal  distribution  occurs  when  the  pressure  is  constant.  As 
mentioned  before,  no  previous  classical  theories  could  predict  how  the  induced  velocity 
would  behave  for  skewed  flow.  Figure  4.33  shows  that  the  induced  velocity  has  a  gradient 
when  in  forward  flight  conditions.  The  rotor  forward  position  is  for  rcos(p  =  —1 ,  where  the 
inflow  distribution  is  smaller  and  it  steadily  increases  towards  the  aft  rotor  position. 


Since  the  general  pattern  for  different  skew  angles  is  similar,  let  us  discuss  the  cases  for 
X=60“.  This  case  is  chosen  because  the  figure  of  merit  has  a  value  half-way  between  the 
ideal  of  1  and  the  minimum  for  axial  flow  ()C=0°).  Starting  with  the  simplest  case  shown  in 
Figure  4.39,  the  rotor  presents  a  constant  pressure  distribution  and  a  sloped  inflow 
distribution,  as  discussed  above.  When  the  combination  of  tilt  and  advance  increases  a 
region  of  downwash  appears  at  the  center  of  the  rotor.  The  downwash  keeps  increasing  and 
shifting  sideways.  This  shift  is  mainly  due  to  the  fact  that  for  skewed  flow,  the  sine  terms  of 
the  equations  are  present  and  contribute  to  the  shape  of  the  pressure  and  inflow 
distributions.  The  sine  part  of  the  definition  for  both  pressure  and  induced  velocity  is  now 
present,  and  it  changes  sign  for  (|j>180°,  where  the  rotor  changes  from  advancing  to 
retreating  side  and  the  reverse  flow  region  occurs. 
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Similar  conclusions  can  be  drawn  for  the  cases  where  x=30°  or  any  other  skew  angle  ranging 
between  these  values.  The  main  differences  are  the  point  at  which  the  downwash  appears 
and  the  magnitude  of  this  downwash. 

Figures  4.51  and  4.52  show  the  pressure  and  induced  velocity  distributions  for  edgewise 
flow  (X~90°)  and  all  values  of  tilt  and  advance  ratios.  For  edgewise  flow,  these  distributions 
do  not  change,  as  was  discussed  previously.  Of  importance  is  to  notice  that,  in  particular, 
the  profile  of  the  inflow  has  a  slope  (similar  to  what  happened  for  all  other  skewed  flow 
cases).  It  is  also  important  to  notice  the  appearance  of  spikes  or  peaks  at  (areas  of  increased 
inflow)  at  (|j=90“  and  (|j=270°.  For  pure  edgewise  flow,  all  the  incoming  flow  is  at  the  plane 
of  the  rotor.  This  means  that  no  flow  will  develop  on  the  perpendicular  plane.  It  is,  in  some 
ways,  similar  to  what  happens  to  a  fixed-wing,  and  the  points  where  the  peaks  occur  are 
really  tip  vortices  due  to  the  fact  that  the  rotor  does  not  have  an  infinite  length.  (Recall  that 
the  flow  is  coming  from  (|j=180°). 

Figures  4.53  through  4.58  show  specific  views  for  all  of  the  cases  in  the  previous  figures 
(including  the  ones  for  axial  flow  shown  in  Section  4.2.3).  Figures  4.53  and  4.56  show  the 
top  view  of  the  rotor  area.  These  are  important  to  show  because  in  the  previous  figures  the 
downwash  appears  to  be  located  in  the  center  of  the  rotor,  while  in  fact  it  is  shifting  towards 
the  reversed-flow  region  as  the  climb  and  advance  increase  (increasing  n).  For  an 
explanation  for  the  reversed-flow  region  see  Appendix  F.  Figures  4.54  and  4.57  show  side 
views  for  the  pressure  and  induced  velocity  distributions  for  the  same  cases.  The  flow  is 
coming  into  the  rotor  at  r  cos  (p  =  —\,  and  r  cos  ^  =  1  is  the  aft  of  the  rotor.  The  asymmetry 
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is  depicted  in  these  figures,  showing  a  downwash  on  the  retreating  side.  This  asymmetry  is 
also  present  on  the  aft  view  shown  in  Figures  4.55  and  4.58,  were  the  downwash  is  clearly 
present  on  the  left-hand  portion  of  the  rotor. 


Pressure  distribution  for  X“30.  u  =  0 


r*cos(ur) 

Figure  4.31:  Pressure  Distribution, =  30,  u  =  0.0,  FM=1.0. 
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Pressure  distribution  for  x=30,  u  =  0.2 
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Figure  4.32:  Pressure  Distribution,  /  =  30,  u  =  0.2,  FM=0.89. 
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Figure  4.33:  Pressure  Distribution,  x  =  30,  u  =  0.4,  FM=0.72. 
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Pressure  distribution  for  u  =  0.6 
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Figure  4.34:  Pressure  Distribution,  x  =  30,  u  =  0.6,  FM=0.56. 
Pressure  distribution  for  x~30.  u  =  0-8 
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Figure  4.35:  Pressure  Distribution,  x  =  30,  u  =  0.8,  FM=0.44. 
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Pressure  distribution  for  X“30.  u  =  1.0 
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Figure  4.36:  Pressure  Distribution,  ^  =  30,  u  =  1.0,  FM=0.36. 


NTInduced  velocity  distribution  for  X~30.  u  =  0 
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Figure  4.37:  Induced  Velocity  Distribution,  y  =  30,  u  =  0.0,  FM=1.0. 
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Induced  velocity  distribution  for  ^  =  0-2 
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Figure  4.38:  Induced  Velocity  Distribution,  =  30,  u  =  0.2,  FM=0.89. 
Induced  velocity  distribution  for  tJ  =  0.4 
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Figure  4.39:  Induced  Velocity  Distribution,  =  30,  u  =  0.4,  FM=0.72. 
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Induced  velocity  distribution  for  X~30.  =  0-6 


r*cos(v|;) 

Figure  4.40:  Induced  Velocity  Distribution, =  30,  u  =  0.6,  FM=0.56. 
Induced  velocity  distribution  for  x~30.  ^  =  0.8 
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Figure  4.41:  Induced  Velocity  Distribution,  i  =  30,  u  =  0.8,  FM=0.44. 
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Induced  velocity  distribution  for  u  =  1  0 
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Figure  4.42:  Induced  Velocity  Distribution, =  30,  u  =  1.0,  FM=0.36. 
Pressure  distribution  for  u  =  0 
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Figure  4.43:  Pressure  Distribution,  ^  —  60,  u  =  0.0,  FM=1.0. 
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Pressure  distribution  for  u  =  0.2 
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Figure  4.44:  Pressure  Distribution,  =  60,  u  =  0.2,  FM=0.94. 


Pressure  distribution  for  x=60,  u  =  0.4 
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Figure  4.45:  Pressure  Distribution,  i  =  60,  u  =  0.4,  FM=0.84. 
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Figure  4.47:  Pressure  Distribution, =  60,  u  =  0.8,  FM=0.62. 


Pressure  distribution  for  X“60,  u  =  0.6 
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Figure  4.46:  Pressure  Distribution,  i  =  60,  u  =  0.6,  FM=0.72. 
Pressure  distribution  for  x=60,  u  =  0.8 
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Pressure  distribution  for  x=60,  u  =  1.0 


r*cosfwl 

Figure  4.48:  Pressure  Distribution,  x  =  60,  u  =  1.0,  FM=0.54. 


V*lnduced  Velocity  distribution  for  x~60,  u  =  0.0 
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Figure  4.49:  Induced  Velocity  Distribution,  y  =  60,  u  =  0.0,  FM=1.0. 
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Induced  velocity  distribution  for  x=60,  u  =  0.2 
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Figure  4.50:  Induced  Velocity  Distribution,  y  =  60,  u  =  0.2,  FM=0.94. 
Induced  velocity  distribution  for  x~60,  u  =  0.4 
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Figure  4.51:  Induced  Velocity  Distribution,  y  =  60,  u  =  0.4,  FM=0.84. 
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Induced  velocity  distribution  for  x^BO,  u  =  0.6 
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Figure  4.53:  Induced  Velocity  Distribution,  y  =  60,  u  =  0.8,  FM=0.62. 
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Figure  4.52:  Induced  Velocity  Distribution,  y  =  60,  u  =  0.6,  FM=0.72. 
Induced  velocity  distribution  for  x=60,  u  =  0.8 


1-0.5 


124 


Induced  velocity  distribution  for  X“60,  u=1.0 
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Figure  4.54:  Induced  Velocity  Distribution,  x  =  60,  u  =  1.0,  FM=0.54. 


Pressure  distribution  for  u=all 
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Figure  4.55:  Pressure  Distribution,  x  =  90,  aU  values  of  u,  FM=1.0. 
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u’Induced  velocity  distribution  for  u=all 
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Figure  4.56:  Induced  Velocity  Distribution,  =  90,  all  values  of  u,  FM=1.0. 
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Pressure  distribution  for  x=60,  u  =  0.0 
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Pressure  distribution  for  x=60,  u  =  0.8 
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Pressure  distribution  for  x=S0,  u  =  0.2 
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Pressure  distribution  for  x=60,  u  =  0.6 
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Pressure  distribution  for  x^^O,  u  =  1.0 
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Figure  4.57:  XY-plane  View  of  Pressure  Distribution,  i  =  60,  various  u. 
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Pressure  distribution  for  x=60,  u  =  0.0 
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Pressure  distribution  for  x=60,  u  =  0.4 
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Pressure  distribution  for  x^^O,  u  =  0.8 
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Pressure  distribution  for  x=60,  u  =  0.2 
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Pressure  distribution  for  x=60,  u  =  0.6 
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Pressure  distribution  for  x=60,  u  =  1.0 
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Figure  4.58:  XZ-plane  View  of  Pressure  Distribution,  =  60,  various  u. 


128 


Pressure  distribution  for  x=60,  u  =  0.0 
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Pressure  distribution  for  x~60,  u  =  0.8 
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Pressure  distribution  for  x=60,  u  =  0.4 
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Pressure  distribution  for  x=60,  u  =  0.2 
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Pressure  distribution  for  X“60,  u  =  1.0 
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Figure  4.59:  YZ-plane  View  of  Pressure  Distribution,  ^  —  60,  various  u. 
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Induced  velocity  distribution  for  X“60,  u  =  0.4 
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Induced  velocity  distribution  for  x=60,  u  =  0.8  Induced  velocity  distribution  for  x=60,  u=1.0 
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Figure  4.60: 


XY-plane  View  of  Induced  Velocity  Distribution,  y  =  60,  various  u. 
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V^lnduced  velocity  distribution  for  %=60,  u  =  0.0 
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Induced  velocity  distribution  for  x=60,  u  =  0.4 
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Induced  velocity  distribution  for  x=60,  u  =  0.8 
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Induced  velocity  distribution  for  x=60,  u  =  0.2 
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Induced  velocity  distribution  for  x=60,  u  =  0.6 
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Induced  velocity  distribution  for  x=60,  u  =  1 .0 
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Figure  4.61:  XZ-plane  View  of  Induced  Velocity  Distribution,  yr  =  60,  various  u. 
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V*lnducecl  velocity  distribution  for  x=60,  u  =  0.0 
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Figure  4.62:  YZ-plane  View  of  Induced  Velocity  Distribution,  x  =  60,  various  u. 
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With  the  results  in  forward  flight,  the  study  of  minimum  induced  power  will  be  complete 


for  any  flight  condition. 


4.3.2.  Finite  Number  of  Blades 


The  theoretical  development  of  the  cases  for  a  finite  number  of  blades  in  skewed  flow 
will  be  shown  in  this  section.  Theory  for  the  cases  of  an  actuator  disk  and  tilted  lift  with  a 
finite  number  of  blades  in  forward  flight  can  be  derived  using  a  modified  system  shown  in 
the  set  of  differential  equations  (3.9)  and  (3.10)  described  in  Section  3.1. 

[K:]{s:}+v[Ej'{a:}=^{C}  (4.73) 

+  e[EJ'  (A'}  =  (4.74) 

The  influence  coefficient  matrices,  ]  the  matrices  defined  by  Equation 

(3.12)  where  the  row  partitions  where  m  Q  have  been  divided  by  two.  The  multiplication 
enables  the  computation  to  become  simpler  due  to  the  fact  that  the  matrices  can  then  be 
divided  into  symmetric  and  asymmetric  parts.  Both  symmetric  and  asymmetric  parts  of  the 
matrices  are  used.  Note  that  the  symmetric  part  provides  only  even  terms  only,  and  the 
asymmetric  part  gives  the  odd  for  y  =  »  ±  1  terms  only. 
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and  are  the  apparent  mass  matrices  also  provided  in  Section  3.1  where 

again,  the  row  partitions  for  m  Q  have  been  divided  by  two,  and  so  on  for  the  rest  of  the 
terms  in  Equations  (4.73)  and  (4.74). 

So  far,  this  formulation  is  valid  for  a  system  with  an  infinite  or  a  finite  number  of  blades  in 
both  axial  and  skewed  flow  conditions  and  for  an  actuator  disk,  or  for  the  case  where  the  tilt 
on  the  lift  is  considered.  The  changes  made  for  each  of  the  cases  has  been  discussed  in 
previous  sections.  Certain  variations  will  be  included  for  this  particular  case  of  skewed  flow 
and  a  finite  number  of  blades. 

First  of  all,  in  previous  cases  where  an  infinite  number  of  blades  was  considered,  the 
system  described  by  Equations  (4.73)  and  (4.74)  was  reduced  to  a  steady  case.  For  the  case 
of  finite  number  of  blades,  this  is  no  longer  a  steady  system.  Moreover,  for  previous  cases  in 
axial  flow,  the  sine-related  equations  all  reduced  to  zero,  but  that  is  no  longer  the  case,  since 
both  cosine  and  sine  contribute  for  the  case  of  skewed  flow  and  finite  number  of  blades. 

Because  of  the  time  dependence  of  the  system,  the  pressure  coefficients,  and 

become  unsteady.  As  a  direct  consequence  of  the  unsteadiness  of  the  pressure,  the  induced 
velocity  coefficients,  and  also  become  unsteady.  Equations  (4.75)  and  (4.76) 

show  the  unsteady  pressure  and  inflow  coefficient  expressions. 

c}={c}.+  Z  [\t:}  co.pT+\c} 

p=Q,2Q,iQ,...  ' 

C')={C).+  Z  ({C'},cos^7  +  {r:')  sin^7 

p=Q,2Q,iQ,...  '  ^  ^  • 
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The  solution  for  the  set  of  differential  equations  is  given  by  the  induced  velocity  coefficients 
as 

<  }  =  {C  }„  +  Z  ({<  cos#  +  {< }  smpJ 

p=Q,2Q,iQ,-  ' 

P=Q,2Q,-iQ,-  ' 

where  Q  is  the  number  of  blades,  m  is  the  spatial  harmonic  number  of  distributions,  p  is  the 
time-wise  harmonic  of  unsteadiness,  a  are  cosine  of  space,  (5  terms  are  sine  of  space,  cp 
are  cosine  of  time  and  sp  are  sine  of  time. 


The  solution  for  the  space  variables  in  time  can  be  obtained  from  the  general  form  the 
set  of  differential  equations  by  the  use  of  harmonic  balance.  The  expressions  are  shown  in 
Equations  (4.77)  and  (4.78). 

\  =  -  V  \ 

[d:J' 

where  [S,-]  =  [E]"'  +^[k:][E][k:] 

2^apjk:][d:j'  [srj' 

where  [6;"]  =  [E]"‘  +^[e;"][E][e;"]  ■ 
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If  instead  of  an  actuator  disk,  the  case  considered  includes  the  tilt  on  the  lift,  then  the 


apparent  mass  matrices  and  must  be  replaced  by  the  corresponding  ones 

using  the  swirl  correction  specified  in  Section  2.6. 

The  optimization  for  this  case  would  be  to  find  the  pressure  coefficients, 
|r”}o  and|r"|  for  both  cosine  and  sine  that  will  provide  the  Goldstein  circulation 

distribution  shown  in  Figures  (4.17)  and  (4.18)  and  translate  these  to  obtain  the  figure  of 
merit.  The  pressure  coefficients  can  be  obtained  from  the  optimum  circulation  distribution 
given  by  reference  [10]  that  uses  a  generalized  dynamic  wake  model  with  a  correction  for  the 
swirl.  Since  it  has  been  proven  that  finite-state  methods  are  able  to  provide  the  Goldstein 
distribution,  the  same  theory  can  be  used  for  forward  flight  with  the  changes  needed  to 
accommodate  for  the  skewed  flow  conditions. 


4.4  Special  Case  of  Hover 


The  previous  results  show  that  Finite-State  methods  can  be  used  to  compute  induced 
power  for  a  rotor  in  cUmb.  It  has  been  shown  that  a  formal  optimization  with  finite-state 
models  was  successfully  performed  in  closed  form  and  that  such  an  optimization  recovers 
the  form  of  solution  for  all  of  the  classical  theories  for  a  lifting  rotor  under  Ughtly-loaded 
conditions.  That  is,  it  agrees  with  the  Prandtl  solution  for  an  actuator  disk  with  a  finite 
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number  of  blades,  with  the  Betz  distribution  for  a  lifting  rotor  with  an  infinite  number  of 
blades,  and  with  the  Goldstein  solution  for  a  lifting  rotor  with  a  finite  number  of  blades. 

According  to  Betz,  the  optimum  induced  flow  distribution  is  one  that  maintains  a 
constant  helix  screw  angle,  (j) .  The  induced  optimum  flow  and  circulation  are  nearly 
identical  to  what  can  be  obtained  from  momentum  theory  and  the  Prandtl  tip  correction 
term,  k.  The  theory  of  Betz  also  applies  to  the  case  of  hover.  The  power  and  thrust  can  be 
determined  from  the  lift  and  drag  perpendicular  and  parallel  to  the  vortex  sheet,  and  can  be 
expressed  in  terms  of  the  lift  and  drag  coefficients  respectively.  For  the  optimum  rotor,  one 
should  find  the  optimum  circulation  that  will  result  in  the  maximum  lift-to-drag  ratio.  If  the 
effect  of  drag  in  the  vertical  direction  is  neglected,  then  one  can  formulate  some  simple, 
closed-form  results  for  the  Uft,  power,  and  thrust  in  hover  of  an  optimum  rotor. 


4.4.1.  Closed-Form  Results  in  Hover 

First  of  all  is  to  define  the  geometry  of  the  problem.  According  to  Betz,  the  optimum 
induced  flow  distribution  is  one  that  maintains  a  constant  helix  screw  angle,  (j).  Thus,  for 
optimum  performance: 


W  —  Wq  cos^ 

(4.79) 

Qx 

cos  (p  —  , 

-y^Q^x^+(V+W()) 

(4.80) 

•  A  V+Wq 

sin  (p  =  1 

-Jq^x^+(V+Wq) 

(4.81) 
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WqQx 

-Jq^x^  +  (V+Wq)^ 


(4.82) 


where  Wq  is  a  nominal  value  chosen  to  obtain  a  given  thrust  for  a  given  rotor  radius  and 
rotational  speed  The  induced  flow  normal  to  the  disk  w^  and  the  swirl  induced  flow  Wg  are 
given  by: 


Wn  =  Wq  cos^  (j)  - 


r^2  2 

WqL2  X 


Q^x^+(V+W() ) 


Wg  =  Wq  cos  ^  sin  (f)  — 


Wq  (V+Wq)Q 

Q^x^  +  (V+Wq)^ 


(4.83) 


(4.84) 


Induced  Flow  computation 

Due  to  the  Betz  geometry  of  similar  triangles  (and  the  fact  that  w  is  always  perpendicular 
to  the  vortex  sheet,  the  following  are  alternative  representations  to  the  relationships  in 
Equations  (4.80)  and  (4.81)  of  cosine  and  sine  of  the  sheet  angle. 

V  cos  (j) 


sin^  = 


cos^zi  = 


^(Qx  -  cos  ^ sin  (f)"-  +  (k  +  cos" 


Qx  +  cos  (j)  sin  (j) 


^(Qx  -  a^Q  cos  sin  +  (K  +  a^Q  cos" 


(4.85) 


(4.86) 


In  some  situations,  the  above  are  more  convenient. 

The  induced  flow  and  circulation  can  be  determined  from  momentum  theory,  the  lift,  and 
the  Prandd  tip  correction  term,  k. 


QL=Qpr 


^  Qx  —  Wq  cos  ^  sin  ^  ^ 
cos^zi 


=  47ikxpwQCOs^^  V+Wq  cos'  (4.87) 


Equation  (4.87)  implies  that. 
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^  ATvkxovf f,  cos 6  ( 2  ,\ 
L  = - V+WoCosV) 


(4.88) 


^  _  ATrkxpwQ  cos^  (j)  V+WqCos^^ 

Q  Q.X  —  W(j  cos  (j)  sin  (j) 

After  substitution  of  the  geometrical  formulas  for  cos  (f)  and  sin  (f) ,  one  obtains 

^  _  Ankp^Q  V+WqQ^x^ 

Q  X  +(V+Wq)  +(V+Wq)^ 

^  _  Ank^Q  V+WqQ^x^ 

’QV+(V+Wof  ^^V  +  V(V+Wo) 


For  lighdy-loaded  rotors  (i.e.,  tPQ«V^,  the  above  formulas  revert  to  the 
relationships. 

Ankp{Q.x^  Vwq  cos (j) 

ne 

4;^/feVW(j  cos^  (j) 

■  Q5 

For  the  case  of  hover,  l/=0,  the  relationships  become: 


_  4;^k/7(Qx)wQCOs^^ 


4;^k/7wQ  Q^'x”* 


Anhpvf^Q  Q^x^ 
Q.Q  Q^x^  +  Wq 


These,  then,  form  the  basis  of  an  efficient  propeller  design. 


(4.89) 

(4.90) 

(4.91) 

Goldstein 

(4.92) 

(4.93) 

(4.94) 

(4.95) 
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Power  and  Thrust 


The  power  and  thrust  can  be  determined  from  the  lift  and  drag  perpendicular  to  and 
parallel  to  the  vortex  sheet.  The  lift  and  drag  on  each  blade  can  be  expressed  in  terms  of  the 
lift  coefficient  and  drag  coefficient. 


L  =  ^  (Qx  — cos^zJsin^)  +{v  +  (j)^  C 


(4.96) 


D  =  (Qx  — cos^zJsin^)  (j)^  C 


(4.97) 


The  thrust  and  power  arise  from  the  geometry  of  the  vortex  angle. 


T  =^[_Lcos^-Dsin^] 


(4.98) 


P  =  Q  (Ox)  [L  sin  (/>  +  D  cos 


(4.99) 


We  assume  that  the  angle  of  attack  at  each  radial  station  has  been  chosen  to  give  the  best 
CJ Cq  and  that  the  chord  has  then  been  chosen  to  give  the  ideal  inflow. 


ir^+v: 


(4.100) 


where  r  is  the  radial  position. 


First,  the  integral  of  the  thmst  results  in  a  relationship  between  and  inflow  ratio,  Vq. 


C,.=2vlB^  l-2Voln  1  +  ^ 


(4.101) 


where 
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B  =  l- 


2ln{2)v, 


(4.102) 


From  the  above,  one  can  find  (for  any  given  C^)  the  appropriate  Vq  for  that  desired  thrust. 
Note  that  we  have  included  the  tip-loss  correction  in  terms  of  Vq. 

Similarly,  the  total  power  coefficient  (including  profile  power)  can  be  done  in  closed  form 

C 

based  on  the  equation  for  the  power,  which  shows  that  the  profile  power  must  be  — —cot(p 
times  the  induced  power.  Thus,  we  have  (after  integration  in  closed  form)  for  hover: 


Cp=2vlB^  l-2vQln  1  +  ^ 


(i  +  Vq^)J  3-^  C 


2  2  2-10Vq"-15Vq"  ^  ,  f  1 

3  Q  1  +  Vq^  ^Kq 


(4.103) 


Note  that,  because  this  is  an  optimum  rotor,  Vq  is  controlled  by  chord.  Thus,  as  Vq  goes  to 
zero,  the  chord  and  profile  power  also  go  to  zero.  The  figure  of  merit  can  be  obtained  from 


the  above  equations  as  — ^ — .  For  induced  power  only,  this  reduces  to: 

V  p 


\  ]  V 

Hover  F.M.  =  B  l-2vQln  1  +  —  + — ^ 

I  ^oJ  1  +  ^0 


(4.104) 


This  is  very  close  to  (but  not  identical  to)  the  climb  figure  of  merit. 


Climb  F.M.  =  BH1-;L"  In  1  + 


(4.105) 


where  the  tip  loss  factor,  B,  is  expressed  in  terms  of  the  non-dimensional  cMmb  rate,  X. 
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Equations  for  and  Cp  imply  an  optimum  chord  (i.e.,  optimum  solidity)  for  maximum 
Figure  of  Merit.  For  small  Vq,  the  profile  power  dominates;  and  Figure  of  Merit  goes  to 
zero.  For  large  Vq,  the  log  term  causes  the  Figure  of  Merit  to  decay  as  1/Vq^.  Thus,  there  is 
an  optimum  thrust,  Vq,  chord,  and  solidity  for  any  Uft-to-drag  ratio. 

As  stated  previously,  for  the  optimum  rotor,  one  should  find  the  combination  of  angle  of 
attack  or  blowing  that  will  result  in  the  maximum  CJ C^.  The  expression  for  optimum  drag- 
to-lift  ratio  is: 


C, 


21n2^  +  4vQMn 

£ 


^1^ 


-2yi 


l-6vQMn 


v^oy 


+  6Vn 


(4.106) 


The  use  of  this  expression  for  a  particular  airfoil  will  yield  a  maximum  drag-to-Uft  ratio. 


4.4.2.  Example  Calculation 


Results  were  obtained  for  an  optimum  propeller  in  hover.  The  results  are  shown  for  the 
case  of  infinite  number  of  blades  and  for  a  propeller  with  four  blades.  Of  importance  is  the 
plot  of  figure  of  merit  variation  with  thrust  coefficient  shown.  The  figure  of  merit  decreases 
for  higher  thrust  coefficient,  and  the  effect  of  finite  number  of  blades  can  be  seen  in  Figure 
4.59. 
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Figure  4.59:  Figure  of  merit  variation  with  thrust  coefficient  in  hover 
(no  induced  drag). 


Figure  4.59  shows  the  case  where  the  effect  of  Cj-,  has  been  neglected.  Figure  4.60  shows  the 
same  graph  for  a  specific  NACA  0012  airfoil  with  a  Uft-to-drag  ratio  of  22,  for  which  the 
effect  of  the  drag  coefficient  is  considered.  It  is  seen  that  there  is  an  optimum  value  for  the 
thmst. 
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NACA  0012  airfoil;  L7D  =  22 


Figure  4.60:  Figure  of  merit  variation  with  thrust  coefficient  for  a  NACA  0012  airfoil. 


The  optimum  Figure  of  Merit  for  different  drag-to-lift  ratios  is  shown  in  Figure  4.61. 

As  it  has  been  shown,  different  quantities  can  be  plotted  to  determine  different  aspects  of 
the  optimum  rotor,  such  as  the  optimal  inflow  ratio  for  different  drag-to-Uft  ratios. 
Moreover,  if  a  specific  airfoil  configuration  is  provided,  results  can  be  shown  for  the  specific 
case  in  climb  and  in  hover  conditions. 
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Figure  4.61:  Optimum  figure  of  merit  change  with  drag-to-Uft  ratio  effect. 


The  above  shows  how  to  design  an  optimum  rotor  (including  the  effects  of  lift  hit,  dp  loss 
and  profile  power)  with  a  special  case  in  closed  form  for  hover.  The  best  possible  figure  of 
merit  is  far  bellow  1.0  but  is  sdU  larger  than  is  being  obtained  today. 

Betz  geometry  was  used  to  obtain  power  and  thmst  for  a  helicopter  rotor  and  applied  to 
the  special  case  of  hover.  The  result  is  a  closed-form  of  the  figure  of  merit  that  can  be  used 
to  determine  how  the  performance  is  affected  by  different  conditions  (finite  number  of 
blades,  lift-to-drag  ratio,  etc.). 
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Chapter  5 


Summary,  Conclusions  and  Future  Work 

5.1  Summary  and  Conclusions 


This  dissertation  shows  various  configurations  can  be  studied  for  induced  power 
calculations  using  Finite-state  methods.  The  development  of  the  solution  is  based  on  the 
need  to  achieve  optimum  induced  power,  therefore  the  optimization  is  set  up  to  achieve  this. 
However,  the  way  in  which  the  optimization  is  performed  allows  for  different  optimizations 
to  be  studied.  Optimum  pressure  coefficients  are  obtained  for  the  minimum  induced  power 
distribution.  However,  this  could  be  changed  to  optimize  for  different  situations.  For 
example,  if  the  need  to  optimize  arises  for  a  particular  blade  shape  with  concrete  dimensions, 
the  vector  containing  the  influence  coefficients  can  be  modified  to  include  the  particular 
characteristics  of  the  blade.  Once  the  theory  has  been  derived,  many  other  cases  can  be 
studied  using  the  same  concepts.  In  a  similar  manner,  the  process  can  also  be  altered  for 
conditions  other  than  that  the  thmst  must  remain  a  constant,  providing  different  results. 

The  results  show  enough  information  in  axial  flow  for  a  finite  number  of  blades  to  help 
in  the  determination  of  how  exactly  a  finite  number  of  blades  affects  the  increase  of  induced 
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power.  For  this  reason,  study  of  the  cases  in  forward  flight  (skewed  flow)  for  a  finite 
number  of  blades  is  unnecessary  at  this  point.  The  theory  has  been  developed  for  this 
particular  case,  but  no  results  have  been  computed.  In  the  future,  some  results  for  figure  of 
merit  or  induced  power  coefficient  calculations  could  be  shown  for  the  sake  of 
completeness. 

A  complete  inflow  theory  is  developed  to  describe  the  flow  in  all  regions  around  the 
rotor.  Previous  cases  for  regions  above  and  on  the  rotor  disk  become  particular  solutions  of 
this  general  formulation.  The  theory  is  then  applied  for  the  case  of  a  quadratic  optimization 
of  power  for  a  constant  thmst.  The  optimization  enables  for  the  figure  of  merit,  or 
efficiency  of  the  rotor,  to  be  found.  The  figure  of  merit  is  shown  to  be  considerably  less 
than  the  ideal  (unity)  for  various  cases,  by  the  calculation  of  optimum  pressure  and  induced 
velocity  distributions  at  various  skew  angles,  increasing  tilt  and  advance  ratio.  The  solutions 
validate  the  theory  for  axial  flow  when  compared  to  classical  results.  In  addition,  it  provides 
a  formulation  to  use  for  ideal  and  real  lightly-loaded  rotors  in  forward  flight.  It  is  found  that 
the  figure  of  merit  decreases  with  advance  and  tilt  due  to  the  downwash  on  inflow  present 
primarily  on  and  around  the  reversed-flow  region. 

The  solution  for  the  optimal  pressure  coefficients  is  performed  and  the  values,  along  with 
the  graphical  representation  of  pressure  and  induced  velocity  distributions  show  that  the 
solution  is  trimmed  forward  to  aft  on  the  rotor  but  that  it  is  not  trimmed  sideways. 
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5.2  Future  Work 


5.2.1  Theory 


Once  the  theory  has  been  developed,  many  applications  can  be  derived  from  the 
complete  formulation.  The  most  direct  line  of  work  that  could  be  developed  from  the 
theory  is  to  apply  it  to  obtain  results  to  see  how  the  formulation  can  be  validated  for  the 
region  below  the  disk.  Also  with  the  complete  formulation,  many  more  applications  can  be 
analyzed.  The  complete  theory  can  be  implemented  into  computer  programs  and  applied 
for  rotor-body  interactions,  ground  effect,  noise  reduction,  etc. 

Another  direct  application  of  the  current  theory  would  be  to  apply  the  complete  inflow 
theory  to  in  forward  flight  for  a  finite  number  of  blades  and  see  how  the  number  of  blades 
affect  induced  power  and  efficiency  in  forward  flight. 

Different  computer-based  optimizations  could  be  used  and  introduced  into  a  CFD 
(Computational  Fluid  Dynamics)  environment  to  obtain  results  for  various  configurations. 
CFD  provides  a  more  accurate  de  description  of  the  flow,  but  it  is  time  consuming  and 
cannot  be  implemented  for  real-time  simulations.  There  are  many  other  computer 
applications  and  comprehensive  codes,  such  as  RCAS  and  CAMRAD  (I  and  II),  simulation 
codes  such  as  FLIGHTLAB,  etc.,  that  would  benefit  from  the  incorporation  of  the  provided 
theory  to  apply  to  rotors  with  specific  characteristics. 
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The  theory  could  also  be  included  into  comprehensive  codes  to  apply  to  real  rotors  and 
blades  and  see  how  the  figure  of  merit  and,  most  importantly,  the  induced  power  are 
affected  and  what  changes  would  be  needed  to  improve  the  efficiency  of  blade 
configurations  used  currently. 


5.2.2  Optimization 


Future  work  in  optimization  could  apply  the  methodology  to  various  constraints  that 
improved  stmctural  hardware.  That  work  could  then  determine  how  much  figure  of  merit 
be  improved  by  various  devices. 

Also  on  optimization,  the  optimal  pressure  coefficients  could  be  constrained  to  obtain  the 
trimmed  solution. 
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Appendix  A 


Ellipsoidal  Coordinate  System 


The  ellipsoidal  coordinate  system  (v,  r|,  y/  )  h  used  in  the  dynamic  wake  model  because  it 
is  convenient  to  express  the  boundary  conditions  for  the  problem  in  this  system. 
Additionally,  a  solution  by  the  use  of  separation  of  variables  for  the  Laplace  equation  takes 
the  form  of  Legendre’s  associated  differential  equation.  The  solution  to  the  separated 
equations  is  easily  recognized  and  applied. 


The  ellipsoidal  coordinate  system  is  defined  as 


where 


X  =  —yjl  —  v^  ^Jl  +  rJ^  cosi^ 

(A.1) 

j  =  yll  —  v^  +  sin^ 

(A.2) 

II 

1 

(A.3) 

-1<V<+1 

(A.4) 

0<7J  <co 

(A.5) 

0<i^<27i 

(A.6) 

Figure  A.l  shows  the  coordinate  system  viewed  in  the  plane  with  the  outline  of 
contours  for  some  v  and  r\  values. 
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Figure  A.l:  Ellipsoidal  Coordinate  System,  xz-plane  view. 


Examples  of  the  contours  of  three  dimensional  plots  for  given  values  of  r\  and  v  are 
shown  on  Figures  A. 2  and  A.3.  Figure  A. 2  shows  the  contour  for  a  fixed  r\  of  0.4.  Figure 
A.3  is  the  contour  on  a  three  dimensional  plot  for  a  v  of  0.8. 
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Figure  A.2:  Ellipsoidal  Coordinate  System  Contour,  rj  —  0.4,  all  v,  all  (//^ 
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Figure  A.2:  Ellipsoidal  Coordinate  System  Contour,  v  =  0.4  for  ^  <  0,  v  =  -0.4  for  ^  >  0,  all  rj,  aU  y/ 
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Appendix  B 

Normalized  Associated  Legendre 
Functions 


Laplace’s  equation  is 


V"O  =  0 


Expressed  in  the  ellipsoidal  coordinate  system  (v,  T),  )  becomes 


A 

dv 


(1 


-V 


do 

dv 


(1  +  77") 


ao 

dtj 


[v^+ri^]  do 
+  dy/ 


Equation  (B.3)  shows  the  separation  of  variables  used  to  solve  the  partial 
equation  of  Equation  (B.2) 


0{ri,v,y/)  =  0^{ri)0^{v)0^{y/) 


Substituting  Equation  (B.3)  into  (B.2)  and  separating  the  variables 


d 

dv 


(> 


-V 


'  dv 


+ 


m 


1-E 


-  +  »(»  +  l) 


Oj  =0 


d 

drj 


(i+'/d 


‘>‘^2 

drj 


+ 


m 


l  +  rj' 


-  — »(«  +  l) 


®2=0 


— +  =  0 
dy/ 


where  m  and  n  are  separation  constants. 


(B.l) 

(B.2) 

differential 

(B.3) 

(B.4) 

(B.5) 

(B.6) 
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Equations  (B.4)  and  (B.5)  are  forms  of  the  associated  Legendre  differential  equations,  and 
the  associated  Legendre  functions  of  the  first  and  second  kind  are  solutions  for  Equations 
(B.4)  and  (B.5)  respectively. 


The  associated  Legendre  function  of  the  first  kind  is 


m  K 


X(-i) 


{2n-2k)\ 


k=0 


2”kl{n-k)l{n-2k)\ 


■V 


-2k 


(B.7) 


The  normalized  associated  Legendre  function  of  the  first  kind  defined  by  Equation  (B.8) 
is  used  throughout  this  dissertation  because  it  provides  a  complete  set  of  functions  that  can 
be  used  to  expand  the  pressure  and  induced  flow  velocity  distributions  radially. 


p;(K)  =  (-ir 


(B.8) 


where  p”  = 


I  1  {n-Vm^. 

2«  +  l  (pi  —  m)\ 


The  normalized  associated  Legendre  function  of  the  second  kind,  (Z;/)  ,  is  explained 

more  in  detail  in  Ref.  2  and  3.  Equation  (B.9)  shows  the  relationship  between  the  associated 
and  the  normalized  associated  functions. 


Q:(in) 


a:  (in) 

fi('O) 


(B.9) 
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Appendix  C 

General  Galerkin  Method 


The  Galerkin  method,  also  known  as  the  method  of  weighted  residuals,  is  an 
approximate  method  to  find  the  solution  of  a  differential  equation. 

Given  a  differential  equation,  L(u)=f,  on  a  domain  D,  with  initial  and  boundary  conditions, 
one  can  find  an  approximate  solution,  that  will  be  dependent  on  x  and  /  and  has  an 
expression: 

V 

=  (C-1) 

7=1 

where  the  (j) ^ ’s  are  the  shape  functions,  which  are  a  complete  set  of  predetermined  functions 

that  have  to  be  chosen.  They  are  also  called  trial  functions  and  must  be  a  linearly 
independent  set  of  functions  and  must  satisfy  the  boundary  conditions  and  the  initial 
conditions  if  any  are  present.  The  goal  is  then  to  determine  the  parameters  a^. 

The  approximate  solution  is  then  introduced  into  the  differential  equation.  Because  the 
solution  is  not  exact,  there  is  an  error  in  the  approximation  that  has  to  be  found.  The  error 
is  called  a  residual  and  can  be  seen  in  Equation  (C.2).  It  is  important  to  reduce  the  residual, 

for  this,  the  error  is  minimized  by  multiplying  the  residual  by  some  test  functions,  and 

then  it  is  integrated  over  the  domain,  D. 

R  =  L{,^)-f  (C2) 
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(C.3) 


jj  =  jj 

D  D 

The  type  of  test,  or  also  called  weighting,  functions  depends  on  the  approximation  that  is 
needed.  There  are  different  types  of  approximations  that  depend  on  what  is  satisfied.  For 
all  of  the  cases,  the  weighting  functions  are  set  to  be  equal  to  the  shape  functions. 

=  (C.4) 

Substituting  Equation  (C.l)  into  (C.3)  yields  an  expression  that  provides  the  (^^  coefficients. 

N  r  1 

S  ll“'*r(y)‘®  (C.5) 

7=1  \_  D  J  -D 

where  N  is  the  fewest  number  of  members  needed  for  a  complete  set  of  functions. 
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Appendix  D 

Transformation  Matrices  in  Closed  Form 


The  closed  form  expressions  for  the  transformation  matrices  were  obtained  by  Ref.  7. 

The  definition  of  the  transformation  matrix  J  comes  from  the  solution  for  the  induced 

flow  expansion  coefficients  in  the  general  set  of  differential  equations  that  define  the 
combined  He  theory  (Eqs.  3.9  and  3.10).  The  definition  for  the  transformation  matrix  is 
given  in  general  form  by: 

(D.i) 


1  _ 

let  'F)  =  —PJ  (k)  and  a  closed  form  expression  for  the  matrix  is: 
V 


—  n+j  2r  2  ]-J  -  -  + 

B'  =B'  =(-l)  2  ^^2n  +  \J2  /  +  1  y  V  ‘Z  1 

jn  jn  \  )  \  Tjr  ^  \  -P  Lu  ? 


q=r,r+2,. 


(»  — ^)(»  +  ^  +  l) 


(D.2) 


The  formula  often  used  is  the  one  for  a  special  case  of  the  inverse,  since  the  expression 
for  requires  more  work. 


[Br  T 

(- 

L 

«+ J—2r 


'Xfi  +  1-^2  y  "1“  1 


(»  +  7)(»  +  7  +  2)  [n-j^ -\ 
and  is  identical  to  for  the  special  case  =  PJ  [v) 

Note  that  [hJ  =  [f]  only  when  n  and  j  go  to  infinity. 


(D.3) 
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Appendix  E 
Table  Method 


In  order  to  compute  the  induced  flow  using  the  generalized  dynamic  wake  model  in 
Reference  2,  one  must  deflne  the  number  of  harmonics  chosen  to  expand  the  induced 
velocity.  The  number  of  harmonics  is  chosen  based  on  the  highest  dynamic  frequency  of 
interest  for  any  given  case.  This  number  provides  the  highest  power  or  Y  shown  in  Table 
E.l.  The  tmncation  for  the  number  of  radial  shapes  is  given  also  in  this  table,  and  depends 
on  the  highest  power  of  Y  for  the  problem. 


HIGHEST 

POWER  OF 

Y 

M 

TOTAL 

INFLOW 

STATES 

0 

1 

2 

3 

4 

5 

0 

1 

1 

1 

1 

1 

3 

2 

2 

1 

1 

6 

3 

2 

2 

1 

1 

10 

4 

3 

2 

2 

1 

1 

15 

5 

3 

3 

2 

2 

1 

1 

21 

Table  E.l:  Choice  of  the  Number  of  Spatial  Modes  (n+m=odd) 
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For  example,  if  the  highest  power  of  Y  is  3,  then  we  will  have  three  harmonics  = 
3).  The  zero  harmonic  («?=0),  will  have  two  inflow  states,  the  first  harmonic  {m—\)  will  have 
2  terms  but  four  inflow  states  (two  for  the  cosine  terms  and  two  for  the  sine  terms),  the 
second  harmonic  {m—7)  will  have  1  term  that  corresponds  to  two  inflow  states,  and  so  on. 
The  total  number  of  inflow  states  for  this  example  will  be  10. 

A  more  detailed  explanation  on  the  choice  of  number  of  inflow  states  to  be  used  for  a  study 
is  given  in  References  [2]  and  [3] . 

In  order  to  determine  the  number  of  inflow  states  needed  to  obtain  accurate  results  for 
the  present  research,  the  figure  of  merit  for  different  number  of  harmonics  was  obtained. 
Table  E.2  and  Figures  E.l  and  E.2  show  the  values  of  figure  of  merit  for  the  simplest  case. 
That  is,  an  actuator  disk  (no  filt)  with  an  infinite  number  of  blades  in  axial  flow  0°).  The 
ideal  figure  of  merit  is  1.0.  However  this  number  is  only  achieved  as  the  number  of  terms 
approaches  infinity. 

Table  E.2  shows  that  the  value  of  figure  of  merit  for  a  case  with  ten  harmonics  (highest 
power  of  T)  is  sufficiently  close  to  the  ideal.  Therefore,  it  can  be  assumed  that  the  use  of 
ten  harmonics  will  achieve  accurate  results.  Twenty  harmonics  provides  a  better  figure  of 
merit  result,  but  the  increase  in  computational  time  might  not  be  worth  the  degree  of 
additional  accuracy.  The  number  of  harmonics  used  in  this  research  is  twenty.  This 
statement  contradicts  the  previous  assumption  of  accuracy  to  be  enough  with  ten  harmonics, 
but  since  this  is  the  first  attempt  to  obtain  results  for  skewed  flow  using  finite-state  methods, 
results  were  obtained  with  the  greatest  accuracy  possible.  However,  in  future  calculations 
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FIGURE  OF  MERIT 


HIGUEST  POWER  OF 


Appendix  F 

Influence  Coefficients  as  a  Fit  of  Coscp 


The  influence  coefficient  vector,  |C"| ,  are  a  Legendre-function  representation  of  cos(p 
which  is  the  loss  of  thrust  due  to  the  tilt  of  the  lift  vector.  Once  more,  they  are  defined  as: 


where 


cl  =  ^ Jo  cos  (v/)  vdvdy/ 

(F.l) 

*=  —  [  [  cos(pP^  {v')vdv  co?,{my/)dy/ 

(F.2) 

C  =— Jo  J  cos^3P„*(K)Wl/sin(/^^^)d'^^ 

n  •'0  •'0 

(F.3) 

cos(p  = 


r  +  ju  sin  Xj/ 


\J{f'  +  jus'mx//y  +/i^ 


(F.4) 


for  a  lighdy-loaded  rotor. 

Then,  the  relationship  between  the  coefficients  and  the  function  coscp  (cosine  of  the  tilt 
angle)  becomes: 

cosg)  =  +  cos{m\j/)  +  C”‘  (F.5) 

n  ^  m,n  ^ 
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This  should  be  true  for  aU  skew  angles,  climb  and  advance  conditions,  and  azimuthal 
locations  on  the  rotor.  To  prove  that  this  is  actually  the  case,  the  distributions  are  plotted 
and  compared.  Figure  F.l  shows  the  distribution  for  a  random  case  using  the  expression  for 
the  cosine  of  the  tilt  in  Equation  (F.4).  Figure  F.2  shows  the  same  distribution  for  a  random 
case  but  using  the  right-hand  side  of  Equation  (F.5)  (i.e.,  influence  coefficient  and  Legendre 
functions  expansion).  Both  of  them  have  the  same  basic  shape  and  dimensional  values. 
Flowever,  the  use  of  the  coefficients  does  not  quite  match  the  close-form  distribution,  as  it  is 
shown  in  Figure  F.3. 


v*cos((t>)  formula  for  |i=0.5196  and  X.=0.3 


■  --0.2 


1-0.4 
-0.6 
-0.8 


r*cos(u») 

Figure  F.l:  v  COScp  Distribution  to  be  fitted  by  Legendre-function  representation. 
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v*cos((t))  fit  for  ^.=0,5196  and  X=0.3 


r’cosfi);) 

Figure  F.2:  Legendre-function  Representation  of  vcoscp. 


Figure  F.3  shows  different  views  for  both  a)  costp  formula  and  b)  Legendre-function 
representation.  The  top  (or  XY)  view  shows  agreement,  but  differences  are  more 
perceptible  when  the  side  views  are  considered.  Figures  F.3a.3  and  F.3b.3  in  particular, 
show  the  tendency  of  the  Legendre-function  representation  (F.3b.3)  to  mimic  the  coscp 
distribution  (F.3a.3).  The  main  difference  is  the  fact  that  Figure  (F.3a.3)  has  a  singularity  for 
t^O,  whereas  its  Legendre-function  counterpart  does  not.  The  Legendre-function 
representation  is  trying  to  introduce  the  singularity,  but  the  number  of  terms  used  for  the 
influence  coefficient  vector  definition  is  not  enough  to  obtain  the  effect.  Twenty  harmonics 
(231  inflow  stated  by  the  table  method  in  Appendix  E)  were  used,  which  has  been  proven 
sufficient  to  obtain  induced  power  and  figure  of  merit  calculations  along  this  dissertation. 
Therefore,  although  the  Legendre-function  representation  shows  slight  disagreement  with 
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what  would  be  thrust  lost  due  to  the  swirl  velocity  induced  by  the  tilt  on  the  lift  for  these 
many  terms,  an  addition  of  a  large  amount  of  terms  would  be  necessary  to  attain  a  closer  fit 
but  would  not  provide  considerable  increased  accuracy  for  neither  figure  of  merit  nor 
induced  power  calculations.  The  provided  results  are  for  a  skew  angle,  of  60°  and  u=0.6, 
which  results  on  a  tdt,  k,  of  0.3  and  an  advance  ratio,  p,  of  0.5196.  Other  cases  where 
analyzed  and  the  same  conclusions  could  be  drawn. 

Finally,  Figures  F.4  through  F.6  show  excerpts  from  Figures  F.l  and  F.2  for  specific 
azimuthal  locations  (specific  angle)  around  the  rotor.  The  figures  shown  correspond  to 
(|j=10°  and  190°  for  Figure  F.4,  (|j=60°  and  240°  for  Figure  F.5,  and  (|j=120°  and  300°  for 
Figure  F.6.  The  solid  lines  represent  the  cosine  of  the  hit  equadon  and  the  dotted  lines  are 
the  Legendre-funcdon  fit  representadon  for  all  radial  posidons.  Figure  F.7  is  a  sketch  of  the 
rotor  viewed  from  above.  To  understand  what  is  happening  and  why  the  fit  for  some  of 
these  distribudons  differs  from  the  close-form  of  coscp  we  dnd  the  locadon  of  Figures  F.4, 
F.5  and  F.6  on  Figure  F.7.  Figure  F.4  corresponds  to  the  distribudon  on  the  crossing  solid 
line  on  Figure  F.7  approximately.  The  differences  are  mosdy  present  on  the  vicinity  of  root 
of  the  rotor  (^0)  since  there  is  a  discontinuity  present.  The  rest  of  the  points  across  the 
rotor  are  adequately  sadsfied  by  the  Legendre-funcdon  representadon.  Figure  F.5 
corresponds  to  the  distribution  on  the  crossing  dashed  line  on  Figure  F.7  approximately. 
Differences  are  greater  in  this  figure  because  the  distribudon  crosses  the  area  of  reversed 
flow. 
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The  reversed-flow  area  is  a  region  on  the  rotor  where  the  flow  is  reversed.  This  means 
that  the  total  velocity  of  the  flow  is  negative  (that  is,  the  flow  has  a  tendency  to  move  on  the 
opposite  direction).  The  area  of  reversed  flow  is  created  by  the  direction  of  the  resulting 
total  velocity  (U-j)  at  a  given  point.  The  incoming  velocity  seen  by  the  blades  in  the  rotor 
(Usifp\i)  is  added  to  the  velocity  of  the  individual  blades  (Qx).  For  180°<(|j<360°  the  incoming 
value  of  the  flow  velocity  is  negative.  Therefore,  if  Usm’\)  >  Q,x  for  a  given  location  x 
(0  <  X  <  1),  the  direction  of  the  total  velocity  is  negative  (reversed  flow).  To  see  where 
these  locations  are  in  the  rotor-frame 

Figure  F.6  corresponds  to  the  distribution  on  the  crossing  dotted  line  on  Figure  F.7 
approximately.  This  location  passes  through  the  complete  reversed-flow  region.  This  is  the 
reason  why  the  fit  shows  the  most  discrepancies. 
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r*sin(\(;) 


cos(4>)  formula  for  |j.=0.5196  and  X=0.3 


cos(4>)  fit  for  |j.=0.5196  and  X=0.3 


Figure  F.3:  Top  and  side  views  for:  a)  costp  Formula  and  b)  Influence  Legendre-function 

representation  of  costp. 
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